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PREFACE 


This  report  presents  the  results  of  improvements  and  extension  of 
the  two-dimensional  distinct  element  code,  UDEC,  and  development  of  the 
data  structure  and  skeleton  code  for  a  new  three-dimensional  distinct 
element  program. 

The  work  was  performed  for  the  U.  S.  Army  Waterways  Experiment 
Station  under  contract  DACA39-82-C-0015 .  These  improvements  and  exten¬ 
sions  of  the  code  supplement  the  original  report  "UDEC  -  A  Generalized 
Distinct  Element  Program  for  Modeling  Jointed  Rock,"  written  by  Dr.  P.  A. 
Cundall  in  March  1980  for  the  U.  S.  Army  European  Research  Office  and 
Defense  Nuclear  Agency  under  contract  DAJA  37-)9-C-0543. 

Mr.  J.  Drake  of  the  Waterways  Experiment  Station  initiated  this  pro¬ 
ject  and  the  final  report  was  prepared  after  consultation  with  Mr.  Drake 
and  Mr.  B.  Armstrong,  also  of  the  Waterways  Experiment  Station. 

Commander  and  Director  of  the  Waterways  Experiment  Station  at  the 
time  of  publication  of  this  report  was  COL  Robert  C.  Lee,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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DEVELOPMENT  OF  GEN ERALIZED  2-D  AND  3- D 


DISTINCT  ELEMENT  PROGRAMS  FOR 


MODELING  JOINTED  ROCK 


PART  I:  INTRODUCTION 


Background 


1.  The  Universal  Distinct  Element  Code  (UDEC)*  is  the  latest 
and  most  advanced  numerical  program  available  for  simulating  the  be¬ 
havior  of  discontinuous  geologic  systems  subjected  to  high  and  tran¬ 
sient  loads.  UDEC  provides  in  one  package  all  of  the  capabilities 
that  existed  separately  in  previous  distinct  element  codes.  The  pro¬ 
gram  is  built  around  a  very  powerful  data  structure  and  is  able  to 
handle  simultaneously  the  interaction  of  a  mixture  of  rock  blocks 
that  have  different  types  of  deformabil ity. 

2.  During  the  initial  development  of  UDEC  several  facilities 
were  encompassed  by  the  original  design  but  were  only  implemented  in 
skeleton  form.  Features  such  as  joint  constitutive  behavior,  dynamic 
cracking,  fluid  flow  and  fluid  pressure  effects  were  identified  as 
requiring  supplemental  work  in  order  to  realize  the  full  modeling 
potential  of  the  code.  Also,  some  utilitarian  improvements  were  sug¬ 
gested:  an  improved  capability  for  dealing  with  flying  blocks  for 
impact-type  problems,  automatic  zoning  for  fully-deformable  blocks, 
improved  logic  for  handling  special  cases  of  splitting  such  as  splitting 
through  corners,  and  more  general  specifications  for  boundary  conditions. 

*  P.  A.  Cundall,  "UDEC  -  A  Generalized  Distinct  Element  Program  for 
Modeling  Jointed  Rock,"  Final  Technical  Report,  European  Research 
Office,  U.  S.  Army,  London,  1980. 


3.  In  addition,  it  was  recognized  that  the  next  logical  exten¬ 
sion  of  the  distinct  element  method  would  be  the  development  of  a  three- 
dimensional  version.  The  first  step  in  this  formidable  task  would  be 
the  design  and  testing  of  a  data  structure  and  test-bed  code  which  would 
be  appropriate  for  three-dimensional  analysis. 


Scope  of  Present  Study 


4.  The  purpose  of  the  present  study  was  to  address  the  consider¬ 
ations  arising  from  the  original  development  of  UOEC.  The  first  ob¬ 
jective  was  to  complete  all  the  unfinished  facilities  identified  above. 
This  accomplished,  the  revised  version  of  UDEC  now  has  a  general  appli¬ 
cation  to  the  following  principal  areas  in  jointed  rock  modeling: 

a.  Discontinuous  systems  can  be  modeled  as  assemblages  of 
blocks  or  particles  of  differing  deformabi 1 ity ;  either 
rigid,  simply-deformable  (with  3  degrees  of  freedom)  or 
fully-deformable  (internally  decomposed  automatically  into 
finite  difference  zones). 

b.  Nonlinear  constitutive  models  including  dilatant  and  non- 
dilatant  behavior  can  be  prescribed  for  both  the  intact 
rock  and  the  discrete  joints. 

c.  Blocks  can  break,  repeatedly,  in  accordance  with  a  user- 
supplied  cracking  criterion. 

d.  Fluid  flow  and  fl uid pressure  generation  in  joints  and 
voids  can  occur  with  flow  rate  specified  in  terms  of 
joint  permeability  and  apparent  aperture. 

e.  Directional  loads  can  be  applied  to  individual  blocks 

and  pressures  can  be  prescribed  to  regions  between  blocks. 

f.  Blocks  or  groups  of  blocks  can  be  explictly  defined  by 
the  user  as  flying  blocks  for  impact  problems. 

5.  The  second  objective  of  this  study  was  to  begin  the  develop¬ 
ment  of  a  new  three-dimensional  distinct  element  program.  A  data 
structure  was  developed  which  was  well -suited  for  the  extension  of  the 
method  to  3-D.  A  test-bed  code  was  then  produced  to  evaluate  various 


aspects  of  the  program  such  as  the  logic  defining  the  characteristics 
of  the  block,  the  detection  of  contacts  and  the  sequence  for  processing 
calculations.  This  effort  has  culminated  in  a  workable  but  primitive 
distinct  element  program  for  three-dimensional  analysis. 

6.  This  report  contains  a  description  of  the  improvements  made 
to  UDEC  and  a  discussion  of  the  development  of  the  three-dimensional 
program.  In  addition,  a  revised  user's  manual  for  UOEC  and  a  new 
user's  manual  for  the  test-bed  3-D  code  are  given  as  appendixes  to  this 


PART  II:  IMPROVEMENT  AND  EXTENSION  OF  UDEC 


Work  Items 


7.  In  the  original  report  several  areas  were  identified  which 
required  additional  work  to  realize  the  full  capacity  of  the  two-dimen¬ 
sional  distinct  element  code.  Specifically,  these  work  items  are: 

a.  complete  edge-to-edge  contact  logic  and  install  a  simple 
constitutive  model  for  rock  joints; 

b.  install  fluid  flow  and  fluid  pressure  generation  logic; 

c.  improve  logic  for  dealing  with  flying  blocks,  i.e.,  blocks 
or  groups  of  blocks  not  in  contact  with  other  blocks; 

d.  install  an  automatic  mesh  generator  for  ful ly-deformable 
blocks ; 

e.  design  logic  to  treat  the  case  of  splitting  through  a  cor¬ 
ner  and  allow  re-entrant  splits  (one  line  crosses  a  single 
block  twice); 

f.  install  dynamic  cracking  including  the  redistribution  of 
forces,  stresses  and  displacements  on  splitting,  and 
create  the  framework  for  user  specified  criteria  for  crack 
development;  and 

<j,.  install  more  general  boundary  conditions. 

Modifications  have  been  made  to  UDEC  to  complete  these  facilities.  In 
Appendix  A  a  revised  UOEC  user's  manual  is  given  which  contains  a  des¬ 
cription  of  the  improvements  made  to  the  code  and  a  complete  set  of  in¬ 
put  commands  and  program  guide.  Sample  problems  are  also  given  which 
demonstrate  the  improvements  made  to  UDEC. 


Data  Structure 


8.  The  program  guide,  given  in  Appendix  A,  contains  the  complete 
contents  of  all  the  groups  in  the  data  structure.  Figures  1  through  5, 


4 


reproduced  from  the  original  report,  show  schematically  the  linkage  of 
these  various  groups  and  should  assist  the  user  in  following  through 
the  program  guide.  Figure  1  shows  the  "linked  list"  arrangement  of  the 
main  data  arrays.  Figures  2,  3  and  4  illustrate  the  conventions  for 
pointers  and  links  in  the  block  data,  domain  data  and  contact  data 
arrays,  respecti vely,  and  Figure  5  shows  the  structural  arrangement  of  re¬ 
dundant  memory  groups.  The  program  guide  and  the  figures  will  assist 
the  user  in  making  any  code  modifications. 
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GROUPS  OF  MEMORY. 


PART  III:  APPROACH  TO  THREE-DIMENSIONAL  MODELING 


Introduction 


9.  This  project  Is  concerned  mainly  with  the  planning  of  a  three- 
dimensional  code  based  on  the  distinct  element  method.  It  is  particu¬ 
larly  important  to  design  the  data  structure  in  a  way  that  anticipates 
how  the  data  will  be  used  during  a  typical  simulation  of  the  behavior 

of  a  blocky  assembly.  Each  physical  quantity  should  be  at  hand  when 
needed,  with  the  minimum  overhead  of  searching,  or  redundant  calcula- 
t ions  . 

10.  Even  though  the  objective  of  the  project  was  to  arrive  at  a 
conceptual  framework  for  future  development,  a  working  program,  called 
D3,  was  written.  The  present  deficiencies  in  D3  are  in  the  areas  of 
contact  detection  and  updating  and  block  creation.  However,  some  as¬ 
pects  of  the  program  are  wel 1 -developed:  for  example,  the  data  struc¬ 
ture;  the  physical  equations  of  motion  and  force-displacement  law;  and 
the  determination  of  volumes  and  centroids  for  arbitrary  blocks. 

11.  Throughout  the  program  03,  functions  or  subroutines  are  used 
to  perform  common  vector  operations.  This  simplifies  the  coding  consid¬ 
erably,  at  the  expense  of  some  increase  in  running  time.  All  vector  and 
tensor  equations  in  this  report  are  expressed  in  component  form,  where 
the  subscripts  i,  j  and  k  range  from  1  to  3,  and  the  Einstein  summation 
convention  applies  for  repeated  subscripts. 

Block  Characteristics 


Geometry 


12.  A  three-dimensional  block  is  defined  by  dividing  its  surface 
into  triangular  faces.  Triangles  are  used  instead  of  arbitrary  poly¬ 
gons  for  the  following  two  reasons. 
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a.  A  surface  is  determined  uniquely  by  specifying  three 
points  in  space.  If  four  or  more  are  given,  the  nature 
of  the  surface  is  undefined  and  ambiguous. 

b.  The  data  structure  is  simplified  if  exactly  three  vertices 
are  associated  with  each  face;  three  memory  locations 

can  be  reserved  in  advance.  Similarly,  exactly  three 
pointers  can  be  provided  to  locate  the  three  adjoining 
faces  to  a  given  face. 

13.  There  is  no  loss  of  generality  by  adopting  the  requirement 
that  the  surface  of  a  polyhedron  be  subdivided  into  triangles.  Any 
arbitrary  shape  can  be  devised  by  using  triangles  as  building  blocks, 
including  blocks  with  concave  regions.  At  present,  ’n  03,  the  vertices 
of  each  face  must  be  given  manually,  but  automatic  surface  zoning  should 
be  possible,  using  the  zone  generation  logic  of  UDEC. 

Volume 

14.  The  calculation  for  block  volume  is  based  upon  Gauss's  diver¬ 
gence  theorem,  given  by: 

9p  =  i  /pn.da  (1) 

}x .;  v  'a 

where  p  is  any  scalar,  vector  or  tensor  variable 
v  is  the  enclosed  volume,  and 

n.  is  the  outward  unit  normal  to  an  element  of  surface,  da. 

If  p  is  defined  as  any  vertex  vector  x. ,  equation  (1)  becomes: 

111  a  l/xi  ni  da 

3xi  v  a 


or,  solving  for  v  using  discrete  areas. 


",  a 


(2) 


where  i  's  the  summation  over  all  surface  elements. 


If  area,  a,  is  planar,  x.  n,  Is  constant  over  the  area. 

For  a  triangular  area  defined  by  the  vectors  zj  ;  and  zj  '  (see  Figure  6) 
the  area  calculation  is: 


where  afc  -  eijk  zfa)  z  (b) 
and  the  unit  normal  is 


3  "  \  la|J  *  \  /ak  *  \ 

j  (eijk  the  Permutation  tensor) 


ai 


"i  2a 

Substituting  equations  (3)  and  (4)  in  equation  (2)  produces 


I  Ix  e  z<*>  z0» 

6  ixk  ei jk  zi  z- 


(3) 

(4) 


3 


(5) 


FIQURE  •  VECTORS  DESCRIBING  TRIANGULAR  AREA  • 

The  volume  associated  with  each  face,  when  defined  by  equation  (5),  rep¬ 
resents  the  volume  of  a  tetrahedron  with  a  base  of  area,  a,  and  apex 
at  the  coordinate  axes  origin.  The  block  volume  is  then  found  from  the 
sum  of  the  tetrahedrons.  To  produce  a  positive  tetrahedron  volume  the 
vertices  defining  a  triangular  face  must  be  ordered  counterclockwise 
when  viewed  from  the  axes  origin. 
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Centroid 


15.  The  centroid  of  the  block  is  calculated  by  recognizing  that 
the  centroid  and  volume  of  each  tetrahedron  are  related  to  the  block 
centroid  by: 


r  =  r  (HUM) 
i  rl  V 
v 


(6) 


where  r.  is  the  centroid  vector  for  the  block 

t 

f 

r  is  the  centroid  vector  for  the  Nth  tetrahedron 
1 

v  is  the  block  volume 

v<N>  is  the  volume  of  the  Nth  tetrahedron 

The  centroid  of  each  tetrahedron  is  calculated  directly  from  the  three 
vertex  vectors  (x^J\  and  x^)  that  define  a  block  face.  The 
centroid  lies  along  the  same  vector  as  the  average  of  these  three 
vectors.  By  simple  integration  techniques  it  can  be  shown  that  the  mag¬ 
nitude  of  the  centroid  is  3/4  of  the  average  vector,  so  that  the  tetra¬ 
hedral  centroid  calculation  becomes: 


(N)  =  v(D  +  x(2)  *  J 3) 


1 


-  x 


1 


1 


+  x' 


i 


(7) 


The  block  centroid  is  then  found  by  using  this  equation  in  equation  (6) 
and  sunning  over  all  tetrahedrons  defining  the  block. 

Radii  of  gyration 

16.  This  calculation  is  Incomplete  In  the  present  version  of  D3. 
Only  dynamic  behavior  is  affected  by  the  moments  of  Inertia,  which  are 
now  taken  to  be  equal,  approximately,  to: 
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(8) 


where  r  is  the  average  distance  from  the  centroid  to  vertices  and  m  is  the 
block  mass. 


Physical  Calculations 


Equation  of  motion 


17.  for  each  block,  the  following  equations  are  integrated  twice 
by  central  finite  differences: 


mu.  ^  imu.  =  IF. 

(9) 

(0  i  (0  i  i 

(10) 

where  u. ,  u^  *  components  of  acceleration  and  velocity 

e’j ,  8 ■  s  components  of  angular  acceleration  and  velocity 
■  sum  of  forces  acting  on  block 
*  sum  of  moments  acting  on  block 
m  «  mass  of  block 


»  moments  of  inertia  about  1,  2,  3  axes, 
a  *  damping  coefficient 


Knowing  the  centroid  motion  and  the  current  locations  of  vertices  and 
centroid,  the  velocities  (and  hence  increments  in  displacement)  of  ver 
tices  are  calculated  as  follows: 


i(p> .  ;<•>>  (,<!>.,<!>> 


i  jk 


(ID 


where  (p)  refers  to  a  vertex 

(b)  refers  to  the  centroid  of  the  block 
e.,,  is  the  permutation  tensor 
’Jk  15 
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The  moment  acting  at  the  centroid  owing  to  a  force  acting  at  a 
surface  point  p  Is  given  by: 


M  z  -  /x(p)  .  y(b)\F(P> 

M1  e1jk(X  j  X  i  ,F  k 


Contact  forces 

18.  At  each  contact,  the  relative  velocity  of  the  two  opposing 
points  is  calculated  using  equation  (11)  for  both  points  and  subtracting 


(c)  '(B)  -(A) 

U  T  '  =  U  j  '  -  U '  •  ' 


where  (A)  and  (B)  denote  the  opposing  points  on  blocks  A  and  B. 

The  relative  contact  velocity  ^  is  resolved  into  normal  and  shear 
partitions: 


where  n.  is  the  contact  normal. 

Normal  and  shear  force  increments  are  then  calculated  as  follows: 


iF(n)  .  k(n)  4t 


aF<‘>  ■  k‘s*  at  -e 


. ,,  e.  F.  n  n' 
ijk  ka8  j  a  B 


where  k^  3  normal  contact  stiffness 
k^  »  shear  contact  stiffness 
nQ  =>  previous  contact  normal 
n^  3  current  contact  normal 


The  second  term  on  the  right  hand  side  of  (17)  corrects  the  current 
shear  force  for  rotation  of  the  contact  normal  during  the  previous  time 
step.  The  expression  is  approximate  only,  and  assumes  that  cos(ao)M. 
The  contact  normal  may  rotate  because: 


a.  the  two  blocks  concerned  have  rotated  about  a  comnon  axis;  01 


b.  the  contact  location  on  one  or  both  blocks  has  changed; 
hence  the  contact  normal  may  have  changed. 


Contact  forces  may  now  be  updated: 


(n ) :=  F(n)  +  AF(n) 

(18] 

(s):=  F(s)  ♦  af(s) 

(i9: 

if  if(5} 


c  ♦  uF^  then 


:=  F(^  (c  +  uF(n)) 
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where  c  =  cohesion 

u  =  friction  coefficient 
:=  means  "replaced  by" 

Note  that  F^  is  stored  in  program  D3  as  a  scalar,  because  the  contact 
normal  is  stored  independently.  However,  F^  is  stored  as  a  vector 
with  components  referred  to  the  global  axes. 

19.  After  calculation,  the  contact  forces  are  applied  immediately 
to  the  two  blocks  comprising  the  contact  (in  a  positive  sense  to  block 
B,  and  in  a  negative  sense  to  block  A).  Equation  (12)  Is  used  to  com¬ 
pute  the  moment  to  be  added,  where  x^  is  the  contact  coordinate. 


Contact  Characteristics  and  Detec t i on 


Prescription  for  contact  normals 

20.  The  blocks  in  UDEC  have  rounded  corners  In  order  to  elimin¬ 
ate  the  singularities,  force-jumps  and  "hang-ups"  associated  with  sharp 
corners.  In  three  dimensions  the  same  idea  is  almost  unworkable,  since 
a  single  spherical  cap  cannot  be  fitted  to  a  vertex  because  it  will  not 
be  tangent  to  all  adjoining  faces.  Some  kind  of  variable-radius  curve 
would  have  to  be  fitted  to  the  vertex.  It  would  have  to  be  tangent  not 
only  to  adjoining  faces,  but  also  somehow  merge  smoothly  with  adjoining 
edges,  which  would  also  be  rounded.  Although  such  a  scheme  may  be  feas¬ 
ible  in  principle,  its  use  would  add  a  large  computing  overhead,  partic¬ 
ularly  in  the  case  of  s imply-deformabl e  blocks,  where  the  angles  at  ver¬ 
tices  are  continuously  changing. 

21.  A  scheme  has  been  devised  that  overcomes  the  problems  with 
sharp  corners,  and  even  resolves  the  ambiguities  present  in  UDEC  for 
very  large  block  overlays.  A  "prescription"  or  rule  is  proposed  that 
furnishes  a  unique  direction  of  contact  normal  to  be  associated  with 
each  point  within  a  block.  Because  two  blocks  must  overlap  in  order  to 
establish  contact,  the  contact  point  must  lie  within  both  blocks.  The 
prescription  is  consulted  to  find  the  average  contact  normal  for  the 
blocks'  internal  point.  Certain  conditions  must  be  fulfilled  by  the 
prescription: 

At  the  surface  of  a  block,  the  prescribed  normals 
must  coincide  with  the  real  normals  (with  jumps  at 
vertices  and  edges). 

b.  There  must  be  a  smooth  transition  in  normal  direction 
from  point-to-point  within  the  block. 

c.  The  rate  of  change  of  normal  direction  with  respect  to 
coordinate  should  reduce  as  the  depth  of  penetration 
Increases. 

In  essence,  the  prescription  provides  a  field  of  normal  vectors  for 
every  internal  point  as  illustrated  ir.  the  figure  below. 
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FlOURE  7  :  NORMAL  VECTOR  CONTOURS  FOR  3-0  CONTACT  DETECTION 


22.  Much  of  the  effect  of  UDEC’s  corner  rounding  is  provided  by 
the  new  scheme  because  there  will  be  a  smooth  transition  as  a  contact 
point  moves  around  a  corner.  Furthermore,  there  is  no  need  to  know 
exactly  which  face  is  providing  support  close  to  a  vertex;  the  known 
normal  determines  the  direction  of  sliding  and  the  direction  in  which 
the  normal  force  increment  is  applied. 

23.  The  following  prescription  for  angles  of  contact  normals  is 
only  tentative.  More  experience  with  its  use  in  D3  is  necessary  before 
it  can  be  accepted  as  being  a  reasonable  analog  of  physical  behavior. 

a.  Select  the  vertex  nearest  the  contact  point. 

b.  Determine  the  normal  distance,  d^,  of  the  contact  point 
from  each  adjoining  face,  N. 

c.  Compute  the  average  normal  face  direction,  weighted  accor 
ding  to  1/d (N ) .  If  the  contact  point  lies  exactly  on  one 
face  N  (d(N)  *  0),  then  the  normal  direction  is  that  of 
face  N. 

d.  The  required  normal  is  the  unit  vector  in  the  computed 
di rection. 

The  prescription  fulfills  the  conditions  previously  set  out,  except  that 
there  will  be  a  slight  change  in  normal  angle  for  deep  penetration  when 
the  "nearest  vertex"  changes. 

Types  of  contact 

24.  Although  six  types  of  contact  can  be  identified  physically, 
only  two  are  necessary  for  complete  support  between  two  blocks. 


Types  of  Physical  Contact 

face  -  face 

face  -  edge 

face  -  vertex  j 

}  necessary  for  support 
edge  -  edge  ) 

edge  -  vertex 

vertex  -  vertex 

Each  of  the  six  physical  categories  can  be  constructed  froin  one  or  more 
combinations  of  face-vertex  and  edge-edge.  These  two  latter  categories 
may  be  termed  "logical  contacts,"  which  are  recognized  by  the  detection 
process  and  in  the  formation  of  the  data  structure.  The  physical  behav¬ 
ior  corresponding  to  the  other  categories  can  be  duplicated  by  knowing 
the  appropriate  areas  and  lengths  of  contacts,  in  the  same  way  that  UDEC 
models  the  physical  behavior  of  an  eoge-to-edge  contact  even  though  the 
logical  contacts  are  of  the  corner-to-edge  form.  D3  does  not  contain  this 
logic  in  its  present  state  of  development. 

Contact  detection 

25.  In  any  code  that  models  interaction  between  arbit  ary  blocks 
or  particles  it  is  necessary  to  avoid  exhaustive  searches  for  those  par¬ 
ticles  that  are  touching  because  the  computer  time  for  such  searches 
increases  as  N^,  where  N  is  the  number  of  particles.  Programs  RBM  and 
SDEM  used  a  "box"  classification  scheme.  Cundall  (1980)  discusses  this 
scheme,  and  its  limitations.  UOEC  uses  a  linked-list  scheme  whereby  a 
block's  contact  candidates  are  found  by  local  search  of  its  surrounding 
domains.  However,  the  two-dimensional  data  structure  of  UDEC  has  no 
convenient  three-dimensional  analog,  as  discussed  in  the  next  section. 

26.  D3  uses  a  scheme  for  which  the  search  time  is  proportional 

to  N,  but  which  is  less  efficient  than  UDEC.  D3  maintains  links  between 
blocks  that  are  near  each  other.  A  given  block  can  then  interrogate  this 
group  of  nearby  blocks  in  order  to  detect  potential  contacts.  The  list 
of  nearby  blocks  is  updated  in  the  following  way.  During  an  "update". 
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a  block  interrogates  not  only  its  local  list  of  neighbors,  but  also 
the  lists  of  its  neighbors.  Blocks  that  are  further  than  a  certain 
radius  are  not  added  to  the  list  {or  are  deleted  if  they  are  on  it  al¬ 
ready),  and  blocks  within  the  radius  are  added.  An  "update"  is  only 
performed  on  a  block  after  it  has  moved  by  some  threshold  distance  since 
its  previous  update.  In  this  way,  updating  of  almost-stationary  regions 
is  avoided. 

Data  Structure 

27.  This  section  describes  the  form  and  use  of  the  data  structure 
in  terms  of  the  pointers  and  connecting  links.  The  complete  content 

of  each  data  array  is  set  out  in  Appendix  B.  Program  D3  is  modeled  close¬ 
ly  on  UDEC  as  far  as  structure  and  operating  logic  are  concerned. 

General  considerations 

28.  The  program  UDEC,  which  models  two-dimensional  block  systems, 
maintains  a  data  structure  with  the  same  topological  form  as  the  physical 
assembly.  The  notion  of  representing  blocks  by  circulating  lists  that 
simultaneously  encompass  the  void  spaces  seems  infeasible  in  three 
dimensions.  It  is  possible  to  have  a  stable  assembly  of  three-dimension¬ 
al  blocks  without  having  an  associated  collection  of  isolated  void  spaces, 
or  "domains";  in  some  three-dimensional  assemblies  it  is  possible  to 
journey  from  one  portion  of  the  void  space  to  any  other  without  needing 

to  pass  between  two  blocks  in  contact.  In  two  dimensions,  the  voids 
can  share  the  same  linked  lists  that  serve  to  describe  blocks.  (A  void 
is  traced  by  following  a  counterclockwise  route,  while  blocks  are 
delimited  by  the  same  list,  but  traced  in  a  clockwise  direction.)  This 
convenient  symmetry  is  not  found  in  three  dimensions. 

29.  Program  D3  embodies,  for  three-dimensional  systems,  a  data 
structure  that  ensures  rapid  access  to  data  as  it  is  needed  during  the  * 
calculation  cycle,  but  the  physical  correspondence  of  UDEC's  data  structure 
is  missing.  This  carries  a  penalty  of  more  time-consuming  searches  for 
contacts  and  increased  difficulty  in  representing  fluid  behavior  in  the 


void  spaces.  Figure  8  shows  the  global  lists  that  link  blocks  and  contacts. 
Block  structure 

30.  For  individual  blocks,  the  data  structure  describes  the  block 
geometry  and  also  permits  the  program  to  jump  from  one  face  to  its  neigh¬ 
bors  directly,  and  from  a  face  to  its  bounding  vertices  directly.  Figures  9 
and  10  illustrate  this  scheme.  Triangular  faces,  apart  from  their  physical 
advantages,  noted  earlier,  lead  to  simplified  data  structures  because 
exactly  three  pointers  suffice  to  link  faces  to  neighbors  and  faces  to 
vertices.  The  connectivity  of  faces  and  vertices  is  specified  completely 

by  the  pointers  provided  in  the  data  array  for  faces,  illustrated  in 
Figure  10.  A  knowledge  of  face  and  vertex  connectivity  is  necessary  for 
an  efficient  scheme  to  detect  and  update  contacts  around  a  block.  The 
data  array  for  vertices  contains  only  coordinate  data,  but  each  block 
has  access  to  a  list  of  its  own  vertices  so  that  coordinates  can  be 
updated  as  the  block  moves.  All  coordinates  are  absolute,  as  components 
are  referred  to  the  global  axes. 

Contacts  and  links  between  blocks 

31.  Global  connectivity  of  the  block  system  is  represented  by  a 
series  of  links  between  nearby  blocks.  When  a  block  system  is  created 
Initially,  these  links  are  established  by  exhaustive  search.  However, 
during  operation,  the  program  can  determine  potential  contacts  by  inter¬ 
rogating  just  those  blocks  In  its  Inmedlate  neighborhood.  In  this  way, 
the  computer  time  needed  for  searching  increases  linearly  with  the  number 
of  blocks,  N,  and  not  as  N^.  The  scheme,  however,  is  not  nearly  as  ef¬ 
ficient  as  that  of  UOEC  because  many  more  potential  contacts  need  to  be 
examined  in  03  for  each  block. 

32.  Contacts  come  In  three  forms:  one  is  a  "degenerate"  form, 
and  the  other  two  correspond  to  "real"  contacts .  A  degenerate  contact 

is  a  simple  link  between  nearby  blocks.  The  memory  taken  by  such  a  contact 
is  much  less  than  that  of  a  real  contact,  but  the  pointers  have  the 
same  locations  as  those  in  real  contacts.  This  permits  both  degenerate 
and  real  contacts  to  be  included  in  the  same  scan.  A  code  number  identi- 
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FIGURE  8  :  GLOBAL  BLOCK  AND  CONTACT  LISTS 


TO  OTHER  FACES 
ON  THIS  BLOCK 


TO  OTHER  VERTICES 
ON  THIS  BLOCK 


TO  OTHER  CONTACTS 
ON  THIS  BLOCK 


CONTACT  ARRAY 


FIGURE  9  :  LISTS  ASSOCIATED  WITH  EACH  BLOCK 


F  :  NEIGHBORING  FACE 


V :  VERTEX 


FIGURE  10  :  POINTERS  ASSOCIATED  WITH  EACH  FACE 


fles  each  type.  Pointers  and  lists  associated  with  each  contact  type 
are  Illustrated  in  Figure  11.  The  two  forms  of  real  contact  are: 
vertex-to-face  and  edge-to-edge.  These  two  categories  are  sufficient  to 
capture  all  types  of  physical  contact,  as  explained  previously. 
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ON  BLOCK  1 


FIGURE  11 


POINTERS  AND  LISTS  ASSOCIATED  WITH 
EACH  CONTACT 


PART  IV:  CONCLUSIONS 


Program  UOEC 

33.  The  two-dimensional  program  UDEC  has  been  considerably 
enhanced:  it  can  now  be  used  to  model  a  wide  spectrum  of  problems  rang¬ 
ing  from  continua  to  dlscontlnua;  from  static  to  dynamic;  and  with  or 
without  pore  fluid  Interaction.  The  utility  of  the  canonical*  data 
structure  has  been  confirmed  by  the  comparative  ease  with  which  the  new 
features  were  installed. 


Program  03 

34.  Considerable  thought  has  been  given  to  devising  a  good  data 
structure  and  physical  idealisation  for  representing  three-dimensional 
block  assemblies.  The  result  is  reported  herein;  much  of  the  scheme 
has  also  been  embodied  in  the  test-bed  program  D3.  In  fact  D3  contains 
a  good  deal  more  than  that  required  by  the  contract:  it  Includes  the 
full  equations  of  motion  for  blocks  and  surfaces,  equations  for  inter¬ 
action  of  contacts,  primitive  logic  for  contact  detection  and  updating, 
and  flxed/free  boundary  conditions. 

35.  It  Is  possible  to  run  very  simple  simulations  with  03  as  It 
stands,  but  the  program  Is  still  only  a  skeleton  code.  The  following 
developments  are  suggested,  in  order  of  priority. 

a.  Test  thoroughly  the  prescription  for  contact  normals, 
and,  if  necessary,  propose  modifications. 

b.  Generalize  logic  for  contact  detection  and  updating,  and 
verify  that  it  will  work  under  extreme  conditions. 

c.  Recognize,  and  treat  correctly,  all  six  categories  of 
contact;  install  corresponding  constitutive  models. 

d.  Add  simply-deformable  logic. 


*A  "canonical"  data  structure  is  "a  model  of  data  which  represents 
the  inherent  structure  of  that  data."  Martin,  J.  (1977),  Computer 
Data-Base  Organization,  Prentice-Hall,  Inc. 
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e.  Install  comprehensive  boundary  conditions: 

1.  stress  tensor 

2.  arbitrary  velocity  prescription 

X.  Allow  blocks  to  split,  dynamically  and  statically; 
~  Include  point-to-point  splitting  law  and  Griffith's 
law  for  simply-deformable  blocks. 

g.  Perform  validation  and  simulation  tests. 


APPENDIX  A:  UNIVERSAL  DISTINCT  ELEMENT  CODE  (VERSION  1.2) 

USER'S  MANUAL 


Introduction 


1.  This  manual  describes  the  latest  Improvements  to  the  Univ¬ 
ersal  Distinct  Element  Code  (UDEC)  and  supplements  the  original  report 
"UOEC  -  A  Generalized  Distinct  Element  Program  for  Modeling  Jointed 
Rock"  written  by  Dr.  P.  A.  «_unda  11,  March  1930,  for  the  U.  S.  Army 
(European  Research  Office,  and  Defense  Nuclear  Agency  under  Contract 
DAJA  37-79-C-0548. 

2.  The  improvements  to  UDEC  were  made  in  the  folioring  general 

areas: 

a.  joint  logic 

b.  fluid  flow 

C.  flying  blocks 

d.  automatic  mesh  generator 

e.  general  splitting  logic 

f.  dynamic  cracking  of  blocks 

£.  generalized  boundary  conditions 
Descriptions  of  these  improvements  and  their  applications  ’n  UDEC  are 
given  in  the  next  section. 

3.  The  modifications  to  UOEC  have  led  to  the  development  of 
several  types  of  constitutive  models  for  the  intact  blocks  and  block 
contacts.  Intact  block  behavior  may  be  defined  by  separate  deformation 
and  fracture  laws,  while  either  point-  or  joint-contact  constitutive  models 
may  be  chosen.  The  different  constitutive  behaviors  are  discussed  below 
and  summarized  in  Table  Al. 

4.  This  manual  also  contains  the  revised  user's  input  commands 
for  UDEC  and  an  updated  program  guide.  Input  and  output  files  are 
presented  for  sample  problems  which  illustrate  the  use  of  the  improve¬ 
ments  to  UDEC. 
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Improvements  to  UDEC 


Joint  logic 

5.  UDEC  recognizes  edge-to-edge  contacts  between  blocks  as  joints, 
and  refers  to  a  constitutive  model  that  works  in  terms  of  stresses 
rather  than  forces.  The  joint  logic  is  used  for  those  blocks  or  joints 
that  are  given  constitutive  number  2  or  5  by  the  user  (see  Table  Al). 

In  addition,  some  or  all  of  the  following  properties  for  joints  should 
be  defined: 

sn^  ....  normal  stiffness  (stress/displacement) 
sSj  ....  shear  stiffness  (stress/displacement) 

Cj  ....  cohesion  (stress/displacement) 
f  .  -  friction  coefficient 

Although  the  joint  logic  may  be  set  for  the  whole  block  assembly,  UDEC 
will  still  refer  to  contact  parameters  under  some  circumstances; 
therefore,  these  parameters  should  also  be  defined.  A  joint  reverts 
back  to  being  a  contact  if  it  no  longer  consists  of  planar  block  faces 
in  opposition.  The  point-contact  logic  is  also  used  if  incremental 
normal  deformation  using  the  joint  parameters  would  be  greater  than  that 
using  the  contact  parameters:  i.e.,  if 

sn  .1  .  <  sn 

where  1.  is  the  length  of  the  joint  and  sn  Is  the  contact  normal 

J 

stiffness. 

Fluid  flow 

6.  Flow  may  occur  between  domains  if  a  differential  pressure 
exists  between  the  domains.  Two  types  of  flow  law  are  used,  depenoing 
on  whether  a  contact  or  a  joint  separates  the  domains.  For  a  contact 
the  flow-rate  is 

q  =  pdiffkc 
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where  i$  the  pressure  difference,  and 

kc  Is  a  permeability  constant,  defined  for  contacts,  for  a 
particular  material  number. 

For  a  joint  the  flow  rule  is: 


where  k j  Isa  permeability  constant  for  joints, 
lj  Is  the  joint  length, 
a^  Is  the  apparent  aperture,  defined  as 


a.  =  max(ares,a  -  strn'/sn.) 

J  J 

where  ares  is  the  residual  displacement  (fully  closed), 
is  the  aperture  for  an  open  joint, 
strn'  is  the  effective  normal  stress,  and 


Snj 


is  the  joint  normal  stiffness. 


The  constants  k^,  ares  and  aQ  are  defined  by  the  user  for  a  particular 
material  number. 

7.  In  one  time-step,  tdel ,  the  adjustment  to  pressure,  pdel , 
in  a  domain  is  as  follows: 


pdel  =  Q(bulkw)(tdel)/Ad 

where  Q  is  the  sum  of  flows  into  the  domain, 

bulkw  is  the  bulk  modulus  of  the  fluid, 

A,  is  the  area  of  the  domain, 
a 

For  a  domain  corresponding  to  a  joint, 


(a.  and  1.  defined  previously).  The  quantities  a^  and  1  ^  are  only 
defined  for  a  joint.  If  constitutive  numbers  2  or  5  are  not  set,  the 
domain  corresponding  to  an  edge-to-edge  contact  will  be  assumed  to  have 
an  area  of  Ad(m1n),  which  may  be  set  by  the  user;  otherwise  it  will 
default  a  small  fraction  of  average  block  areas.  For  regular  domains, 
Ad(min)  is  the  limiting  area  for  fluid  calculations. 
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8.  A  printout  of  fluid  flow  in  all  joints  and  contacts  may  be 
requested  by  giving  the  PRINT  FLOW  command.  Aperture  and  length  are 
also  printed  for  joints. 

9.  The  Influence  of  a  fluid  pressure  gradient  is  Included  in 
UDEC  for  fully-saturated  blocks  subjected  to  gravity  loading.  This  has 
been  accomplished  by  adding  a  buoyancy  force  term  to  the  law  of  motion 
for  a  block.  The  buoyancy  force  is  defined  by  a  fluid  density  parameter, 
rhow,  in  the  FLUID  input  command. 

Flying  blocks 

10.  UDEC  keeps  track  of  "flying“  blocks  (i.e.,  blocks  not  in 
contact  with  other  blocks)  by  retaining  one  link  to  the  main  data 
structure.  This  ensures  that  new  contacts  will  be  detected  in  the 
domain  containing  the  flying  block.  The  single  link  is  of  the  same 
form  as  a  regular  contact,  but  it  contributes  no  forces  and  is  deleted 
immediately  after  the  block  comes  into  contact  with  other  blocks. 

Groups  of  flying  blocks  are  handled  in  an  identical  manner.  The  same 
logic  ensures  that  the  group  is  linked  to  the  global  data  structure  by 
one  virtual  contact.  81ocks  which  are  initially  not  in  contact  with 
other  blocks  must  be  linked  to  the  main  data  structure  using  the  LINK 
input  command. 

Automatic  mesh  generator 

11.  The  automatic  mesh  generator  is  based  upon  that  described 
in  the  report  "Computer  Modeling  of  Jointed  Rock  Masses"  written  by 
Dr.  Cundall,  et.  al . ,  (see  Technical  Report  N-78-4  for  the  U.  S.  Army 
Engineers  Waterways  Experiment  Station,  August,  1978). 

12.  Automatic  mesh  generation  for  a  fully-deformable  block  is 
accomplished  in  three  stages.  First,  all  corners  of  the  boundary  are 
linked  so  that  the  block  is  discretized  as  a  triangular  finite- 
difference  mesh.  Then,  the  triangles  are  split  until  all  triangular 
sides  are  smaller  than  a  maximum  edge  length  specified  by  the  user. 
Finally,  all  internal  grid-points  are  adjusted  until  their  coordinates 
coincide  with  the  average  of  the  coordinates  of  the  surrounding  grid- 
points.  The  generator  appears  to  be  sufficient  for  discretizing  most 
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blocks  provided  the  aspect  ratio  (longest  to  shortest  dimensions)  of 
the  block  Is  kept  smaller  than  2:1. 

General  splitting  logic 

13.  The  logic  for  splitting  blocks  has  been  overhauled  so  that 
a  split  may  occur  at  any  orientation.  Splits  through  corners  are 
allowed.  If  a  given  split-line  passes  too  close  to  a  corner,  the  line 
Is  diverted  so  that  the  corner  is  split.  The  criterion  of  "closeness" 
is  based  on  the  given  rounding  length;  the  line  is  diverted  if  a  newly- 
created  corner  would  interfere  with  the  existing  corner  (i.e.,  their 
rounding  arcs  would  overlap).  After  diverting  a  proposed  split-line 
through  corners  (if  necessary)  a  check  is  made  to  see  if  the  line  would 
coincide  with  an  existing  edge;  if  it  would,  the  split  is  rejected  for 
that  block. 

14.  Block  splitting  is  accomplished  via  subroutine  XYSPL( MAT, ICONS). 
This  routine  only  requires  two  coordinates  ((x3,y3)  and  (x4,y4))  which 
define  the  split  line  through  the  block  and  MAT  and  ICONS  which  assign 
joint  properties  and  constitutive  behavior  to  the  newly  created  joint. 
Dynamic  cracking  of  blocks 

15.  UDEC  has  been  modified  to  allow  dynamic  cracking  of  rigid  and 
simply-deformable  blocks.  The  decision  to  check  for  cracking  is  made  by 
introducing  a  tensile  strength  factor,  tf,  to  the  material  property  list. 

If  a  block  has  a  specified  tensile  strength  factor,  it  Is  searched 

once  every  cycle  for  conditions  which  would  satisfy  a  user-supplied 
cracking  criterion.  If  this  criterion  is  met,  the  block  is  split  Into 
two.  The  joint  created  by  splitting  a  block  will  take  the  material  and 
constitutive  numbers  of  the  block. 

16.  Two  cracking  criteria  are  presently  available  in  the  code. 

A  criterion  based  upon  a  relationship  developed  from  "point-load"  testing 
has  been  assigned  to  constitutive  numbers  1  and  2.  The  tensile  strength 
factor  In  this  case  is  defined  by: 

tf  - 

20 

where  fl  and  f2  are  two  opposing  contact  forces  applied  to  the  block,  and 
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d  is  the  distance  between  these  forces.  Splitting  of  the  block  occurs 
If  the  maximum  value  of  the  contact  force-distance  relationship  equals 
or  exceeds  tf.  Dynamic  cracking  is  not  permitted  through  corners  or 
too  close  to  corners  (d  <  1/2  smallest  block  edge)  for  this  cracking 
criterion. 

17.  A  criterion  based  on  Griffith  theory  is  assigned  to  constitu¬ 
tive  numbers  4  and  5.  This  criterion  evaluates  conditions  for  cracking 
In  terms  of  the  Internal  stresses  In  SDEF  blocks  (tensile  stresses  are 
assumed  positive).  The  relationships  for  block  splitting  are  defined  by: 

if  3spl  +  sp2  >  0 


if  3spl  +  sp2  <  0 

where  spl  is  the  maximum  principal  stress  in  the  SDEF  block, 
sp2  is  the  minimum  principal  stress  in  the  SDEF  block, 
tf  corresponds  to  the  uniaxial  tensile  strength  of  the  intact 
material . 

When  stress  conditions  exceed  the  tensile  strength,  the  block  Is  split 
through  its  centroid  in  a  direction  parallel  to  sp2  and  the  block 
Stresses  are  set  to  zero. 

18.  It  should  be  noted  that  these  two  cracking  models  do  not 
account  for  energy  lost  In  the  system  when  the  fracture  occurs.  A  .rore 
thorough  approach  should  take  into  account  the  change  of  train  energy 
into  kinetic  energy  at  failure. 

Generalized  boundary  conditions 

19.  Two  types  of  boundary  conditions  can  be  used  in  UOEC.  X  and  Y 
directional  loads  can  be  added  to  block  centroids  using  the  LOAD  command. 
Domain  pressure  can  be  user-controlled  using  the  PF I X  and  PFREE  commands. 
Summary  of  constitutive  models 

?0.  Each  constitutive  number  gives  the  user  a  different  combina¬ 
tion  of  constitutive  behavior  for  the  intact  block  and  the  contact 
between  blocks.  cour  combinations  are  presently  defined  (see  table  below). 
Other  combinations  are  left  to  the  discretion  of  the  user. 


and 


tf  =  spl 
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Table  A1 

Constitutive  Behavior  Models 


Intact  Block 

Contacts 

Constitutive 

Number 

Deformation 

Law 

Cracking 

Law 

Deformation 

Law 

1 

elastic-isotropic 

point- load 

point  contact 

2 

elastic-isotropic 

point-load 

joint  contact 

3 

elastic-isotropic 

Griffith 

point  contact 

4 

elastic-isotropic 

Griffith 

joint  contact 

Input  Comnands 


Notes:  Upper-case  letters  in  a  command  or  parameter  must  be  typed;  the 
remaining  letters  are  optional.  Lower-case  parameters  stand 
for  numeric  values.  Integers  must  be  given  for  parameters  start¬ 
ing  with  i,  j,  k,  1,  m,  n.  Real  numbers  may  be  given  as  integers, 
but  not  vice  versa. 

Input  is  free-format:  parameters  may  be  separated  by  any  number 
of  the  following  characters,  in  addition  to  spaces: 

»  (  )  .  / 

An  END  command  is  required  at  the  end  of  the  input  file  (after 
the  STOP  command).  The  first  command  must  be  START  or  RESTART. 

*  =  comment  line 

♦  =  continuation  line 


Block  Material  n  Constitutive  m  xl  y 1  x2  y2  ... 

Create  a  rigid  block  of  material  number  n  and  constitutive  number 
m.  Oefaults  are  n=l ,  m=l,  if  m,  n  are  omitted.  Corner  coordinates 
are:  (xl.yl),  (x2,y 2)  etc.,  in  a  clockwise  direction.  Continua¬ 
tion  lines  may  be  used,  but  a  pair  of  numbers  defining  a  corner 
must  not  be  separated.  Only  one  BLOCK  command  may  be  used  per 
run  at  present.  Further  blocks  may  be  created  with  a  SPLIT 
command,  and  unwanted  ones  deleted  with  the  DELETE  command. 

Any  blocks  may  be  changed  to  simply-  or  ful ly-deformable  with 
a  CHANGE  command. 

CHange  xl  x2  yl  y2  Sdef  Material  n  Constitutive  m 

Fdef 

All  blocks  with  centroids  lying  within  the  range  xl<x<x2  ,  yl <y<y2 
are  changed  to  simply-  or  ful ly-deformable  (Sdef  or  Fdef  respec¬ 
tively).  Material  and  constitutive  numbers  may  also  be  changed. 


Oo  n  time-steps  (cycle  0  is  permitted  as  a  check  on  data). 
DAmping  fcrit  freq  Mass 

Stiffness 

Internal 

Viscous  damping  is  applied  in  the  form  of  Rayleigh  damping, 
fcrit  is  the  fraction  of  critical  damping  and  freq  Is  the  center 
frequency.  If  a  qualifier  is  not  given  as  the  third  parameter, 
full  damping  is  used.  The  word  "Mass"  eliminates  the  stiffness- 
proportional  dashpots.  The  word  "Internal"  causes  the  specific 
damping  to  be  applied  to  the  3  internal  degrees  of  freedom  of 
simply-deformable  blocks. 

OElete  xl  x2  yl  y2 

All  blocks  are  deleted  in  the  range  xl<x<x2  ,  yl<y<y2 
Oump  n  m 

Dump  memory  to  printer  from  the  main  array  from  address  n  to 
address  m.  Internal  pointers  MFREE,  JUNK  IBPNT,  ICPNT  and  IDPNT 
are  also  printed.  MFREE  gives  the  highest  memory  location  that 
is  currently  free. 

End 

Last  input  command. 

Fix  xl  x2  yl  y2 

All  blocks  are  fixed  in  the  range  xl<x<x2  ,  yl<y<y2. 

FRAction  f 

f  is  taken  as  the  fraction  of  critical  time-step  to  be  used. 
Fluid  rhow  bulkw 

Fluid  properties  are  specified  for  an  effective  stress  analysis. 
These  are  the  density,  rhow,  and  the  bulk  modulus,  bulkw. 

FRee  xl  x2  yl  y2 

All  blocks  are  set  free  in  the  range  xl<x<x2  ,  yl<y<y2. 

Note:  By  default  all  blocks  are  free  initially. 

Generate  xl  x2  yl  y2  Manual  Gridpoints  (glist)  Zones  (zlist) 

Automatic  (amaxl) 


All  blocks  encountered  In  the  range  xl<x<x2  ,  yl<y<y2  are 
discretized  as  ful ly-deformable.  For  automatic  generation  the 
parameter  (amaxl)  must  be  given  to  define  the  maximum  edge 
length  of  the  triangular  zones.  For  manual  generation  a  list 
of  grid-points,  (glist),  and  zones,  (zlist)  must  be  given. 

The  format  for  (glist)  is: 

xl  yl  x2  y2  x3  y3 . . 

where  each  x,y  pair  Is  a  coordinate  of  a  grid-point.  The  format 
for  (zlist)  is: 

11  ml  nl  12  m2  n2  . 

Each  triple  corresponds  to  the  three  grid-points  that  define 
the  zone,  where  the  numbering  of  the  grid-points  refers  to  the 
order  in  (glist),  starting  with  the  last  point  (i.e.,  the  last 
grid-point  is  number  1).  The  grid-points  should  be  given  in 
clockwise  order  around  the  zone.  Both  (glist)  and  (zlist)  may 
extend  over  an  arbitrary  number  of  continuation  lines,  but 
doubles  and  triples  should  not  be  split  over  two  lines.  If  a 
given  coordinate  lies  within  a  certain  tolerance  of  a  block 
corner,  the  grid-point  is  placed  on  that  corner.  The  tolerance 
is  taken  as  0.9  times  the  rounding  length.  Grid-point  coordinates 
can  be  defined  to  coincide  with  block  corners  but  should  not  be 
defined  to  lie  along  block  edges,  for  manual  generation. 

Gravity  gx  gy 

Gravitational  accelerations  are  set  for  the  x-  and  y-  directions. 
Link  xl  yl  x2  y2 

Links  a  flying  block  to  the  main  data  structure.  (xl,yl)  are 
the  coordinates  of  any  point  inside  the  flying  block  and  (x2,y2) 
are  the  coordinates  of  any  point  Inside  the  block  which  will 
provide  the  link  to  the  flying  block.  This  block  should  be  the 
one  which  Is  topologically  closest  to  the  flying  block. 

(xl,yl'  and  (x2,y2)  should  be  chosen  close  to  the  blocks' 
centroid  locations  to  ensure  correct  linkage. 

Load  xl  x2  yl  y2  xload  yload 

All  blocks  with  centroids  lying  within  the  range  xl<x<x2  ,  yl<y<y2 
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are  prescribed  static  loads  applied  at  the  block  centroid. 

PFIx  la  p 

The  pressure  Is  controlled  In  the  domain  with  address  la.  The 
real  constant  value  for  pressure,  p,  is  Inserted  In  the  pore 
pressure  offset  of  the  domain  list. 

PFRee  la 

The  pressure  Is  not  controlled  In  the  domain  with  address  la. 

PLot  Nofix  Zones  NC  Vel 

If  no  parameter  follows  the  PLot  command,  all  blocks  and 
centroids  are  plotted.  If  "Nofix"  is  used,  no  fixed  blocks  are 
plotted.  The  keyword  "Zones"  is  used  to  plot  the  zones  In  fully- 
deformable  blocks.  The  word  "NC"  deletes  corner  rounding  on  all 
blocks  and  "Vel"  plots  block  velocity  vectors  at  block  centroids. 
Print  Blocks  Contacts  CORners  Domains  List  DLlst  Flows 

Data  are  printed  on  blocks,  contacts,  corners,  domains  and  linked 
lists  for  blocks  and  domains.  Fluid  flows  in  joints  and  contacts 
are  printed  with  the  FLOWS  keyword. 

PROperty  Material  n  keyword  value 
n 

The  first  parameter  must  be  the  specification  of  the  material 
number.  Material  properties  are  defined  for  material  number  n. 
Property  keywords  are: 


Bulk(or  K) 

bulk  modulus 

G 

shear  modulus 

Density 

density 

KN 

contact  normal  stiffness 

KS 

contact  shear  stiffness 

Cohesion 

contact  cohesion 

Friction 

contact  friction  coefficient 

JKN 

joint  normal  stiffness 

JKS 

joint  shear  stiffness 

JCoh 

joint  cohesion 

JFric 

joint  friction  coefficient 

Tf 

tensile  strength  factor 

JPerm 

joint  permeability  constant 

CPerm  contact  permeability  constant 

AZero  aperture  for  tero  normal  stress 

ARes  residual  aperture  at  high  stress 

(Units  of  joint  normal  and  shear  stiffness  and  joint  cohesion 

are  stress/displacement.) 

Restart 

The  program  is  restarted  using  data  from  the  restart  file. 

RSet  v  la  loff 

The  real  value  v  is  Inserted  in  the  main  array  at  address  la, 
with  offset  ioff. 

ROund  d 

Each  block  corner  is  rounded  with  a  circle  that  is  tangential  to 
the  two  corresponding  edges  at  a  distance  d  from  the  corner. 

SAve 

The  current  problem  state  is  saved  on  the  restart  file. 

SPlIt  xl  yl  x2  y2  Material  n  Constitutive  m 

All  blocks  in  the  path  of  a  line  extending  from  point  (xl.yl) 
to  (x2,y2)  are  split  into  two.  The  joint  created  by  the  split 
is  assigned  a  joint  material  number  n  and  a  joint  constitutive 
number  m.  If  MAT  or  CONS  are  omitted,  the  joint  or  contact 
will  take  the  material  and  constitutive  numbers  of  one  of  the 
adjoining  blocks  (however  no  number  will  be  printed  when  the 
PRINT  CONTACTS  command  is  given.) 

STArt 

The  program  does  a  cold  start. 

Stop 

The  run  stops. 

View  1x1  1x2  iyl  1y2 

The  Integer  ranges  1x1  to  1x2  and  Iyl  to  1y2  define  the  viewport 
region  on  the  plotting  device  within  which  the  plot  will  be  made. 
Oefaults  are  ixl=0,  ix2=2000,  1yla0,  1y2*1400. 

Window  xl  x  2  yl  y2 

The  coordinate  ranges  xl  to  x2  and  yl  to  y2  define  In  real  problem 
units  the  region  of  the  model  to  be  plotted.  Defaults  are 
xl=0,  x2=10,  yl=0,  y2=7. 
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Program  Cuid« 


Pirimtttn  and  Data  Croup 


Offset s  for  block  data  array 


Not#:  The  first  integer  in  each  block  array 
-  (offset  0)  is  th«  block  type  number,  as  follows 

1  rigid  block 

2  simp  I y-def ormab I e  block 

3  ful I y-deformable  block 

KB  Pointer  to  nett  block  in  block  list. 

KP  Pointer  to  one  corner  in  block's  corner  list. 

KMAT  Material  number. 

KCONS  Constitutive  number. 

KBCOD  Code  number: 

0  free  block 
1  fixed  block 

KX  x  coordinate  of  centroid. 

KY  y  coordinate  of  centroid. 

KXD  x  velocity. 

KYD  y  velocity. 

KTD  Angular  velocity  < counterc I ockwise  positive). 
KAREA  Block  area. 

KBM  Block  mass. 

KB  I  Moment  of  inertia. 

KBFX  x  centroid  force-sum. 

KBFY  y  centroid  force-sum. 

KBFT  Centroid  moment  sum. 

KXL  x  load  applied  to  block  centroid. 

KYL  y  load  applied  to  block  centroid. 

KBEX  Extension  pointer  (to  SDEF  or  FDEF  data) 

Offsets  for  corner  data  array 


Notei  The  first  integer  (offset  0)  contains 
-  the  value  MCOR  to  denote  a  corner. 

KL  Pointer  to  next  corner  or  contact  on 
block,  in  clockwise  direction. 

KR  Pointer  to  next  corner  in  counterc I ockwise 
direct  ion. 

KNB  Pointer  to  host  block. 

KXP  x  coordinate  of  corner. 

KYP  y  coordinate  of  corner. 

KXCP  x  coordinate  of  local  circle  center. 

KYCP  y  coordinate  of  local  circle  center, 

KRAD  Radius  of  local  circle. 


KXDP 

KYDP 

KCP 


x  velocity 
y  velocity 
Pointer  to 


of  corner* 
of  corner, 
corresponding 


grid-point  if  block 


is  ful  I  y-deformab  le* 
Offsets  for  contact  data  array 
Note  i 


KCAP 

KCQ 

KCL 


The 

the 


first  integer  (offset 
value  MCON  to  denote  a 


0)  contains 
contact. 


KC 
KB1 
KB2 
KL 1 

KL2 

KD1 

KD2 

KCM 

KCC 

KXC 

KYC 

KXDC 

KYDC 

KCS 

KCN 

KCFS 

KCFN 

KCCOD 


Pointer  to 
Address  of 
Address  of 
Pointer  to 
of  block  co 
Same  as  KL1 
Address  of  domain  to 
going  from  block  KB1 
Address  of  domain  to 
going  from  block  KB2 
Material  type  number 
Constitutive  number, 
x  contact  coordinatt 
y  contact  coordinate 
Relative  x  velocity 


next  contact  in  contact  list, 
first  block  involved  in  contact, 
second  block  involved  in  contact, 
next  item  in  clockwise  list 
rrespond i 
,  but  for 


ng 

to  KB  l  . 

bl 

ock  KB2 

>  1  e 

ft  of  c 

to 

KB2. 

>  r  i 

ght  of 

!  to 

KB1  . 

!  • 

(Of 

block 

contact  i 


KB2  relative 


KB1  ) . 
y  velocity. 


shear  displacement, 
normal  displacement. 


(compression  positive). 


to  block 
Relative 
Relative 
Relative 
Shear  force. 

Normal  force 
Code  number; 

1  corner /corner  contact 

2  corner/edge  contact  ( KB1 ...  corner, 

KB2. • . edge ) 

3  edge/corner  contact  (K81»..edgef 

KB2. . . corner  I 

Mean  aperture  for  joint 
Flow-rate  in  joint  or  contact 
Length  associated  with  joint 


Offsets  for  domain  data  array 
Note  t 


The 

the 


f  irst 
value 


integer  (offset  0)  contains 
MDOM  to  denote  a  domain. 


KD 

KDAR 

KPP 

KUMAX 

KDLOOP 

KDCOD 


Pointer  to  next  domain  in  domain  list. 
Domain  area. 

Pore-pressure  for  domain. 

Fictitious  domain  displacement. 

Pointer  to  one  contact  in  counterc 1 ockwise 
list  around  domain. 

Code  number i 

0  domain  pressure  not  controlled 
1  domain  pressure  controlled 
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Simpl y-dtformib )e  extension  array 


KED1 1  ) 

KED12  )  Strain-rate 
KED21  )  tensor 
KED22  > 

KSI11  ) 

KSI12  )  Internal  stress 
KSI21  )  tensor 
KSI22  ) 

KSAll  )  Applied  stress 
KSA12  )  tensor  (multiplied 
KSA21  >  by  block  area) 
KSA22  ) 

Offsets  for  grid-point  data 


KC  Pointer  to  next  grid-point  in  grid-point  list. 

KCOR  Pointer  to  correspond ing  block  corner. 

KXG  x  coordinate. 

KYG  y  coordinate. 

KXDC  x  velocity. 

KYDG  y  velocity. 

KCFX  x  force -sum. 

KCFY  y  force-sum. 

KGPM  grid-point  mass. 

Offsets  for  zone  data 


KZ  Pointer  to  next  zone  in  zone  list. 

KZC  Start  of  triple  pointer  to  3  surrounding 

gr id-points. 

KZSlt  ) 

KZS12  )  Stress  tensor 
KZS22  ) 

KZM  Zone  mass 

KZLL  Pointer  to  neighboring  zone  for 

mi xed -d iscret ization  calculation. 

Logical  unit  numbers 


LUNIF  Unit  number  for  input  file. 

LUNOF  Unit  number  for  output  file. 

LUNG  Unit  number  for  general  I/O  (e.g.  restart). 

LUNP  Unit  number  for  plotted  output. 


Number  of  words  in  data  arrays 


NVCR 

Corner 

NVBL 

Block 

NVCN 

Contact 

NVDO 

Domain 

NVSD 

Simp  1 y-deformab 1 e  extension 

NVZO 

Zone 

NVCP 

Cr id -point 

Arr  ay 

l  imi  ts 

MTOP  Size  of  main  array  <IA). 

IMMAT  Maximum  number  of  materials. 

NCONS  Maximum  constitutive  numbers. 

NT  i'P  Numoar  of  block  types  (rigid.  SDEF.  etc.) 

Head  codes  (contents  of  first  integer  in  data  groups* 


MR  I C  =  1  Rigid  block 

MSDEF  =  2  Simpl y-deformab le  block 

MFDEF  *  3  Full y-deformab I e  block 

MCOR  Corner 

MCON  Contact 

MDOM  Domain 

Main  Common  Block  Variables 


LINE(60> 
LINEl (30) 
LPNT1 I ) 

RAFLAC 

PPFLAC 

ERFLAC 

STFLAC 

DCF LAC 

COFLAC 

NCFLAC 

CRFLAC 

JMPSAV 

NERR 

JUNK 

MFREE 

I  BLOCK 


Buffer  for  current  input  line  in  A1  format. 

Buffer  for  next  input  line. 

Pointer  to  start  of  parameter  I  in  L1NE(  ) 
after  removal  of  blanks,  etc. 

.TRUE,  if  pore -pressure  calculation  requested. 
•TRUE,  if  an  error  has  occured. 

.TRUE,  if  the  first  input  line  has  been  processed . 
•TRUE,  if  the  domain  pressure  is  controlled. 

•TRUE,  if  the  current  line  is  a  continuation. 
•TRUE,  if  the  next  line  is  a  continuation. 

•TRUE,  if  block  splitting  calculation  is  requested 
Index  of  last  computed  GOTO  in  MON. 

Error  number. 

Pointer  to  list  of  spare  memory  groups. 

First  unused  memory  address. 

Current  block  number. 


I  DOM 

I  STACK 

NCYC 

NCTOT 

TDEL 

FRAC 

I ROUTE 

NLINE 

NPACE 

JMPCEN 

ALPHA 

BETA 

CONI 

CON2 

BDT 

ALPB 

C1B 

r->  o 

DECRAD 

PI 

DAM  IN 
ATOL 

BTOL 

CTOL 

DTOL 

DTOL2 

ETOL 

FTOL 

CTOL 
HTOL 
IBPNT 
ICPNT 
IDPNT 
IODPNT 
AKN ( I ) 
AKS ( I ) 
AMU  < I ) 
COH( I) 
AKNJ ( I ) 
AKSJ ( I ) 
AMUJ < I ) 
COHJ ( I ) 
PERMJ ( I ) 
PERMC ( I ) 
AZEROI I ) 
ARES ( I ) 


contact 


at  which  a  contact 


Currant  domain  number. 

Stack  pointer. 

Currently  requested  number  of  cycles. 

Total  number  of  cycles. 

Time -step. 

Requested  fraction  of  critical  time-step. 
Routing  number,  used  in  main  routine. 

Output  line  count. 

Output  page  count. 

Routing  number  for  continuation  line  in  CEN. 
Mass  damping  coefficient. 

Stiffness  damping  coefficient. 

Damping  factor  ( 1 . 0-ALPHA*TDEL/2. 0) 

Damping  factor  (  1 . 0/ (  1 . 0+'ALPHA*TDEL/2. 0  )  ) 
BETA/TDEL 

Internal  mass  damping  coefficient  for 
s i mp I y - de f ormab I e  blocks. 

Damping  factor  (  1 . 0 - ALPB*TDEL/ 2 . 0 ) 

Damp ing  factor  ( 1 . 0 / < 1 . 0 ♦ ALPB* TDEL/ 2 .  0 ) ) 

PI  /  130 
3. 14159 

Minimum  domain  area  allowed. 

Distance  between  particles  at  which  a  contac 
is  first  formed. 

Distance  between  particles  at  which  a  contac 
is  broken. 

Maximum  (negative)  overlap  allowed 
when  forming  contacts. 

Rounding  length. 

DTOL/2.0  (maximum  contact  overlap) 

Limit  on  maximum  domain  displacement 
to  trigger  contact  update. 

Total  area  of  blocks  for  setting 
plotting  scale  factor. 


Pointer  to  list  of  blocks. 

Pointer  to  list  of  contacts. 

Pointer  to  list  of  domains. 

Pointer  to  outer  domain. 

Normal  contact  stiffness,  material  I. 
Shear  contact  stiffness,  material  I. 
Contact  friction  coefficient,  material 
Contact  cohesion,  material  I. 

Joint  normal  stiffness,  material  I. 
Joint  shear  stiffness,  material  1. 

Joint  friiction  coefficient,  material  : 
Joint  cohesion,  materia)  I. 

Joint  permeability  constant,  material 
Contact  permeability  constant,  material 
Initial  aperture,  material  I. 

Residua)  aperture,  materia)  I. 


a  I  I  owed 


DAM  IN 

Minimum  domain 

area  tor 

f  lui 

DENS( I ) 

Density#  mater 

i  a  I  1 . 

BULK!  I ) 

Bulk  modulus# 

material 

I. 

SHEAR ( 1 1 

Shear  modulus# 

materia  1 

I . 

TFAC( I  I 

Tensile  streng 

th  factor 

#  mat 

ALAMl < 1 1 

Lame  constant# 

material 

I. 

ALAM2 ( I ) 

Lame  constant# 

materia  1 

I . 

CRAVX 

x  component  of 

grav i tat 

iona  1 

CRAVY 

y  component  of 

grav i tat 

iona  1 

RHOU 

Fluid  density. 

BULKW 

1X1 

1X2 

iyi 

IT2 
RX 1 
RX2 
RY  1 
RY2 
I A  (  ) 


Fluid  bulk  modulus. 

Plotter  viewport  coordinate. 
Plotter  viewport  coordinate. 
Plotter  viewport  coordinate. 
Plotter  viewport  coordinate. 
Prob I em  window  coordinate. 
Problem  window  coordinate. 
Problem  window  coordinate. 
Problem  window  coordinate. 


Main  array. 

Main  Subroutine  Calling  Map 


UDEC 

-SETUP 

-MON 

-HALT 

-PRINT 

-CREATE 

-SPLIT 

-APLOT 

-INI 

-CYCLE 

-PPSCAN 

-BLKSCN 

-PPCEN 

-PPDIS 

-CONSCN 

-CRKSCN 

-DOMSCN 

-CEN 
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d  ca I cu I  a t ions. 

erial  1 . 

acce I erat ion. 
acce I erat ion. 


Sample  Problems 


The  following  four  sample  problems  Illustrate  the  Improvements 
made  to  UOEC. 

No.  1  Single  point-load  cracking 
No.  2  Pressurized  cavity 
No.  3  Complex  block  deformation 
No.  4  Projectile  breaking  beam 

The  printed  output  for  each  problem  should  be  used  to  provide  a  check  that  the 
program  is  performing  correctly. 


ft 

o 


START 

FRCP  *T*1  D£W*2000  W-=lM  KS*:E4  F=2  TF=i.5e4 

WAV  0  -10 

Dtff  .1  15  (HASS) 

FRACO.l 

BLOCK  (0,0)  <0,30/  (40,30)  «0,0) 

SPLIT  -1,10  41,10 
C^IT  -1,20  41,20 
SPUT  -1,5  4i,5 
SPLIT  25,2.5  25,12.5 
DELETE  25,40  5,10 
SPLIT  30,15  30,40 
DELETE  30,40  2C.3C- 
*lt®  0  50  j  30 
FIX  0,4/  0,5 
r-LOT 

crc  400 

PftlHT  BLOCKS  COT4CTS 

PLOT 

STOP 

END 


1START 


tmm  UDEC:  IMIVBSAL  DISTINCT  ELBBff  CODE  (VESICH  1.2) 


THIS  IS  A  STMT  AM 

m  MM  D06>2000  W1E8  I3ME4  F*2  TF*2.5£4 
NMW«  -10 
>Dwr  .1  IS  (WES) 
miTNESS-ONf  DC  TOW  SCT  TO  ZERO 

>mc  o.i 

>IUK1  (0,0)  (0,30)  (40,30)  (40,0) 

>9PllT  -1,10  41,10 
ISPIIT  -1,20  41,20 
WFUT-1,5  01,5 
)Sn.IT  25,2.5  25,12.5 
>DELETE  25,40  5,10 
ISPIIT  30,15  30,40 
)DOm  30,40  20,30 
)HM  0  50  0  50 
>nxo,4o  0,5 
)pujt 
)Crc  ooo 

INITIAL  TWESTO>  >  1.000E-02 
CURRENT  CTCU  COM  •  400 

>Ht[)ff  BLOCKS  Off  ACTS 


SUCH  DATA 

HOC*  NAT  CONST  CENTROID  COORDS.  NASS  POL.  NON,  X,T,WETA  VELOCITIES  X,T,DCTA  FORCES 


151 

(RIGID) 

1 

1 

3.4446*01 

1. 140E*01 

2.5946 *05 

4.0558*04 

-1.2138-04  -2.747E-04 

1.904E-05  4.0026*02 

2.5776*04 

2. 4446*03 

1 

(RIGID) 

1 

1 

1.5086*01 

2.4731*01 

4.000£*05 

5.000E+07 

2.9926*05  5.592E-05 

5.4O3E-07  7.784E*02 

4.0017*04 

I.214E.04 

433 

(RICID) 

1 

1 

1 .2506*01 

7.4271*00 

2.500E*05 

1.354E*07 

4.0396-05  *1.3486-04 

4.7096-07  -4.573E.02 

2.500E*O6 

2*7526*03 

314 

(RICID) 

1 

l 

2.0006*01 

2.5006*00 

4.000E*05 

5.417E*07 

0.000E«O0  O.OOOE^OO 

0.0006*00  4.7876*01 

•1.44VE+07 

4.545E*07 

74 

(RIGID) 

1 

1 

1.3556*01 

1.5006*01 

5.3916*05 

3.777E*07 

3.3176-0*  -1 .0026-05 

3.5876-04  -V,OOOE*CO 

5.3876*04 

•7 ,4806*02 

OMTCT  DATA 

CONTACT  NAT  COST 

X,Y  COORDS. 

FORCE  DISftACDWT 

NORM. 

sew  KRW. 

SHEAR 

2*5 

0 

0 

3.8886*01 

4.9956*00 

1.8846*44 

-3.1798*45  -9.8446-03 

3.1796-03 

144 

0 

0 

2.9041*01 

1.9486*01 

3.0416*04 

•4.290E+05  -1.7536-01 

4,2906-03 

41V 

0 

0 

2.4876*01 

9.7596*00 

7.1901*05 

-2.9158*45  -4.9546-05 

1.1246-02 

552 

0 

0 

2.4016*01 

9.8706*00 

4.7156*04 

-1.8171*45  *0.1741 -42 

1.8171-03 

382 

0 

0 

2.4746*01 

1.0176*01 

0.0006*00 

0.0006*40  1. 5446-01 

0.0006*40 

588 

0 

0 

2.4506*01 

4.9586*00 

8.4846*04 

1.5576*45  -8.4846-02 

-1.5576-03 

404 

0 

0 

5.0056-01 

4.9496*00 

4.1156*04 

1.5406*45  -4.1156-02 

-1.5406-03 

243 

0 

0 

5.4326-01 

1.9886*01 

2.9396*44 

3.7406*45  -2.9136-02 

-3.7406-03 

144 

0 

0 

5.0426-01 

9.9146*00 

4.4746*44 

1.4216*45  -4.4246-02 

-1.4216-03 

>fW 

(STOP 
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Sample  Problem  No.  2 


Upper  block  Is  forced  Into  a  cavity  by  an  applied  load.  Pressure 
Is  thereby  Induced  In  the  cavity,  driving  the  rlghthand  block  outwards. 
The  pressure  also  Induces  flows  In  the  surrounding  joints,  and  hence 
pressure-drops  In  the  enclosed  volumes  between  blocks.  The  outer  domain 
Is  held  to  a  fixed  pressure  of  zero. 


START 

PROP  MM  DENS1 2000  M.QE6  CU.OEfl  KX-1.E06  KSM.E08  F*0.5 

PROP  MM  CPEWME-9 

FLUID  0.0  1.019 

DM*  .5  16.  MSS 

FMC  0.10 

MACK  0.,0.  0. ,30.  60. ,30.  40., 0. 

HU®  0.2 

SPLIT  -l.,tO.  41., 10.  MM  C0K3-1 
SPLIT  lb. ,9.  15., 31.  MM  CCNS=1 
SPLIT  25.  ,9.  25., 31.  MT*1  OHM 
SPLIT  10,20  41 ,20  MM,  COhSU 
DELETE  15., 15.  10., 20. 

FIX  0.,40.  0.  ,10. 

FIX  0. ,15.  0. ,30. 

FIX  25. ,40.  20. ,30. 

PFIX  69  0.0 

LOAD  15. ,25.  20. ,30  0.0  -1.0E6 
WIND  0  50  0  40 
CYCLE  3000 

PRINT  FLOWS,  DOMIHS,  BLOCK'; 

PLOT 

STOP 

EHD 


PROBLEM  NO.  2  INPUT  FILE 
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. n.  »rt  «  mU  a.i't  «i ^  | 

WSSmSSmSS* 


gg M 


PROGRAM  UKC>  UNIVERSAL  DISTINCT  ELMXT  COM  (VERSION  1.2) 


THIS  IS  A  START  AM 

HKGP  MM  006*2000  M.0E8  C*1.0E0  BM.E08  XSM.E08  F>0.5 
HOT  MM  CPBIME-9 
)FUIU>  0.0  1*0C9 
>BMT  .5  10.  MOB 

STIFFICSS-DAMIMC  TSOI  SET  TO  ZERO 
>FMC  0.10 

)UOCI  0 .,0.  0.^0.  40. /SO.  40. ,0. 

HUM  0.2 

NPUT  -l.,10.  41.  ,10.  MM  C0NS*1 
)9FUT  15,  ,9,  13.  ,31.  MT*1  <»8«i 
)8FUT  25., 9.  25.  ,31.  MT*1  C0NS*1 
>SPUT  10,20  41,20  MM,  C0NS*1 
IDELETE  15. ,25.  10. ,20. 

>FIX  0.,40.  0.,10. 

)PIX  0.  ,15.  0.^0. 

>FIX  25., 40.  20. ,30. 

>FTIX  49  0.0 

HOAD  15., 25.  20. ,30  0.0  -1.0E6 

XtIM)  0  50  0  40 
> CYCli  3000 

INITIAL  TDCSTEP  <  8.944E-03 
CURRENT  CYCLE  COJNT  *  3000 
MM  FLOWS,  DOWNS,  BLOCKS 

FLOW  ACROSS  CONTACTS  OR  JOINTS  ... 


CONTACT 

X 

Y  FLOW 

LENGTH 

APERTURE 

DOM 

DOC 

M3 

1.500E*01 

1.9I5E*01  Z.944E-05 

0.0006*00 

O.OOOE*00 

522 

905 

704 

2.5006*01 

2.0206*01  2.944E-05 

0.0006*00 

0.0006*00 

522 

728 

440 

3,9006*01 

2.000E*01  -2.921E-05 

0.0006*00 

0.0006*00 

49 

470 

424 

2.4Z5E*01 

2.0006*01  2.9446-05 

0.0006*00 

0.0006*00 

522 

470 

500 

2.425E+01 

1.0006*01  *2.9446*05 

0.0006*00 

0.0006*00 

100 

522 

449 

2.5006*01 

2.0756*01  -2.921E-05 

0.0006*00 

0.0006*00 

49 

728 

350 

1.4006*01 

1.000E*01  2.9446*05 

0.0006*00 

0.0006*00 

522 

372 

205 

1.5006*01 

2.075E*01  -2.921E-05 

0.0006*00 

0.0006*00 

49 

905 

144 

3.9006*01 

1.0006*01  -2.921E-05 

0.0006*00 

0,0006*00 

49 

188 

144 

2.0006*01 

1.0006*01  -2.921E-05 

0. 0006*00 

0.0006*00 

49 

372 

DOMAIN  DATA 

DOMAIN 

PORE  PMSSUDE 

VOLIJK 

HAXMH  DISPLACEMENT 

905 

2.9226*04 

4.800E-01 

0.512E-04 

728 

2.9226*04 

4.0006*01 

5.9515*03 

470 

2.9226*04 

4.0006*01 

3.4506-04 

522 

5,844**04 

9.9996*01 

1.1945*02 

372 

2.9226*04 

4.0006-01 

0.0006*00 

188 

2.9226*04 

4.000E-01 

3.4506*04 

49 

(OUTER 

0.0006*00 

4.8C0E-01 

M40E-02 

BOUNDARY) 


noa  MM 


noa  wt  cost  oflwiD  cotws.  mss  m.  nan. 

558  I  1  3 . 250€<-01  2.500E.01  3,0006*05  8.125 {*08 

(rich) 

311  1  1  3.3551*01  1.5006*01  3, OOOC .05  «.  1251*04 

(RICH) 

185  1  1  7.500E*CX>  2.000E*01  8.0006*05  3.1251*07 

(RICH) 

78  1  1  2.0006*01  2.3V5E«OI  2.000€*05  3.J33£*08 

(RICH) 

1  1  1  2.0006*01  5.000E*00  8. OCX* *05  l.iJ3E*08 

(RICH) 

)K0T 

>STOP 


X,T, THETA  VELOCITIES  X,r,TWTA  FCHS 

O.OOOE*CO  0.0006*00  0.00(*«00  8.5086*03  - 3.2X1*03  1, 5551*03 

3.8571-02  -2.5111-08  -7.4281-11  -I. 1711-01  -2.284-02  -5.4321-02 

0.0008*00  0.000£*00  0.0006*00  -3. 7 141*04  -1.8586*04  -1.7086*05 

3.2186-05  -3.8456-02  7.5586-08  3.0856*04  2.1836*04  J.312E«05 

0.0006*00  0.0006*00  0.0006*00  5.2886-02  -1.0571-01  -14206*00 
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Sample  Problem  No.  3 


A  small  heavy  block  sits  on  a  large  block  that  has  low  moduli 
and  Is  fully  deformable.  After  1000  time-steps  the  plot  shows  the 
complex  deformation  pattern  that  develops,  and  the  printout  gives  the 
Internal  stresses. 


a.  After  1,000  cycles 


Figure  A3.  UDEC  Sample  Problem  No.  3 


START 

PROP  MM  0*2000  JN*1E0  KS*lE0  FRIC«.i  C*iE4  HJLK=2£4 

PROP  HAT *2  0*10000  W*1E0  IS*1E0  FRIC*.l  G*1E6  RJLME8 

ROUND  *  0.2 

BLOCK  2,2  2,0  4,0  6,2 

SPLIT  0,7  7,7 

SPLIT  4.5  4  4.5  9 

DELETE  2  5  7  0 

SPLIT  5.5  6  5.5  9 

DELETE  5.5  6  7  0 

SPLIT  0.3  8,3 

:^ANCE  4.5  5.5  7  6  HAf=2 

'‘NAfttt  5  a  3  7  FDEr 

GEh  2  *  s'  4oTQ: i * i 

Fix  2  >  23 

GFAv  0  -10 

DWf  .1  i  HASS 

CfCLE  1000 

xlriD  0  6  0  10 

PLOT 

FftlKT  &L0-1KS 

STCF 

EHD 


PROBLEM  NO.  3  INPUT  FILE 
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>ST*T 


WOtWl  UCECi  UNIVERSAL  DISTIHCT  ELEMENT  COM  (VBtSION  1.2) 


TVlIS  IS  A  START  HR 

>WflP  MTM  0*2000  DRIES  KS*IES  »IC*.l  C<IES  UU>2E6 

>pwp  mr>2  e-ioooo  dries  ks*ies  fric*.i  c*ies  suu-ies 
hsu#  *  0.2 

>sun  2,2  2,S  4,4  4,2 
JSPtIT  0,7  7,7 
> SPLIT  4.5  4  4.5  9 
)  DELETE  2  5  7  6 
/SPLIT  5.5  4  5.5  9 
>D£LfTE  5.5  4  7  S 
>SPl!T  0,3  4,3 
>OMNCE  4.5  5.5  7  S  HAT«2 
•OWNCS  2  4  3  7  ntt 
>C£K  2  4  3  7  AWOM.t 
/FIX  2  4  2  3 
•CRAV  4  -10 
>0#f  .1  1  !W3S 

ST;ffNESs-:*wPi*  teka  set  to  zero 
•CiCIE  1000 

INITIAL  TIHESTEF  *  1.76vE-03 
CUtfENT  CTCLE  CONT  *  1000 
•MIND  0  4  0  10 
•PLOT 

•PRINT  SLOCKS 


SLOCK  DATA 

SLOCK  (WT  CONST  CEXTkOID  COWS.  NASS  POL,  It*.  X,T,!>€TA  VELOCITIES  X,T , THETA  fOKES 

195  1  1  4.000£*00  2.5001*00  3.000£*03  1.133£*04  0.0006*00  0.000£*00  0.000£*00  1.624E*04  -J.772E*05  -2.10Z£«05 

IRICID) 

76  2  1  5.3771*00  7.374E*00  1.000£*04  1.647E.03  1.I57E-4I  4.<O0E-02  I.01SE-0I  -3.3l#E*03  3.307E9O4  -l.)54E*04 

(RICID) 

I  1  1  4.000E«00  5.000E+00  3.200E*04  S.533E«04  0.0006*00  O.OOOEtOO  O.O0OE*OO  O.OOOE*00  O.OOOE*00  O.OOOE*00 

(FDEF) 


ZONE 

CPID-FOINTS 

S1C71A- 11  SIOW- 12  SIOW- 22 

MASS 

1120 

463 

705 

574 

9.261E*03  2.734E*03  4.3372*03 

5.0006*02 

1129 

429 

574 

705 

1.29SE*04  1.364E*03  -4.0106*03 

5.0006*02 

1136 

463 

442 

705 

2.307E«03  6.611E*03  -6.164E«03 

5.0006*02 

1147 

429 

705 

442 

5.936E*03  5.457E«03  - 1.6666*04 

5.0006*02 

555 

447 

750 

574 

1.7551*04  3.954E«02  -4.475E*04 

5.0006*02 

1154 

463 

574 

750 

-5.3991*03  3.390E*O3  -3.Z32£*04 

5.0001*02 

1145 

447 

415 

750 

1.504£*04  1.1971*03  -3.917E*04 

5.0006*02 

1174 

463 

750 

415 

-7.935E*03  4.1946*03  -2.477E*4M 

5.0006*02 

723 

429 

771 

574 

1.4591*04  -7.903£*03  -2.75SE*03 

5.0006*02 

1163 

525 

574 

771 

2.044E*04  -3.4291*03  -3.4536*04 

5.000E*02 

lift 

429 

433 

771 

1.094E*04  -|,499E*04  -2.9S6E+04 

5.C00E*02 

1201 

525 

771 

433 

1.  >761*04  -1.0701*04  -6.204E*O4 

5.0006*02 

732 

525 

626 

576 

6.329E*03  -4.4741*03  -6.45lE*04 

5.0006*02 

1210 

447 

574 

626 

2.155E*03  -2.497E*03  -5.CHJE*04 

5.0006*02 

1219 

525 

451 

626 

1.0641*04  -  3.722£*03  -7.159E*04 

5.000E*02 

1226 

447 

626 

651 

4.5271*03  -1.724E*03  -5.605E*O4 

5.0006*02 

741 

463 

697 

>42 

7.773E*03  1.1346*03  -5.551E*03 

5.000C*02 
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ZONE 

GRtt-POINTS 

sicnA'ii 

SIGM-12  S1UM-22 

mss 

1237 

411 

442 

697 

1.736E+03 

1.4446*02  7.7106*03 

5.0006*02 

1244 

463 

495 

697 

1.2226*04 

3.2256*03  -3.3026*02 

5.0006*02 

1253 

411 

897 

495 

8.2146*03 

2.2376*03  1.2926*04 

5.0006*02 

404 

454 

607 

615 

4.3036*03 

-1.473E+03  -1.638E+04 

5.0006*02 

1244 

463 

415 

607 

1 >9271*04 

2.0576*03  -1.543E+04 

5.0006*02 

1273 

454 

495 

607 

-9.764E+C3 

2.2606*03  -2.0446*04 

5.0006*02 

1262 

463 

607 

495 

5.3916*03 

5.6066*03  -1.7796*04 

5.000E*02 

760 

447 

424 

615 

1.5456*04 

9.1676*03  -3.9276*04 

5.0006*02 

1291 

544 

415 

424 

2.0456*04 

1.0576*04  -3.7246*04 

5.0006*02 

1300 

447 

449 

424 

4.390E+03 

1.7706*04  -4.9406*04 

5.0006*02 

1309 

544 

424 

449 

9.3466*03 

1.9096*04  -4.7456*04 

5.0006*02 

769 

544 

945 

615 

1.2436*04 

4.2746*03  -4.003E+04 

5.0006*02 

1316 

456 

415 

945 

•1.2096*03 

6.3956*03  *3.397E*04 

3.0006*02 

1327 

544 

494 

945 

4.7236*03 

7.4446*03  -2.6506*04 

5.0006*02 

1334 

454 

945 

694 

•9.3236*03 

9.5826*03  -2.2436*04 

5.000E*02 

796 

402 

924 

433 

5.7486*03 

4.647E*03  -4.4636*04 

5.000E+02 

1345 

525 

433 

924 

1.291E*04 

2.7846*03  -4.5876*04 

5.0006*02 

1354 

402 

436 

924 

5.75l£*03 

5.4496*03  -4.0746*04 

5.000E+02 

1343 

525 

924 

436 

1.3006*04 

1.5276*03  -4.1526*04 

;.0OoK*O2 

837 

445 

656 

436 

8.1996*03 

-1.8076*03  -7.794E+04 

5.000E*02 

1372 

525 

438 

656 

1.1036*04 

-4.421E+03  -4.724E+04 

5.000E+02 

1381 

443 

651 

856 

1.242E*04 

-1.4286*03  -8.098E*04 

5.000E+02 

1390 

525 

856 

451 

1.5746*04 

-4,406E*03  -4.978E+04 

5.0006*02 

647 

445 

666 

651 

1.3576*04 

3<447E*03  -7.943E*04 

5.000E*02 

1399 

565 

651 

668 

7.3646*03 

2.1056*03  -8.0866*04 

5.0006*02 

1406 

445 

467 

666 

8.441E+03 

1.2336*04  -1.1726*05 

5.0006*02 

1417 

565 

666 

467 

1.4046*03 

1.0776*04  -1.1926*05 

5.0006*02 

904 

565 

1031 

651 

•2.144E*04 

1.078E*04  -9.2856*04 

5.0006*02 

1426 

447 

651 

1031 

1 . 129E*03 

3.5856*02  -4.442E+04 

5.0006*02 

1435 

565 

449 

1031 

-3.458E*03 

1.518E+04  -8.915£*')4 

5.000E+02 

1444 

447 

1031 

44 9 

1.904E*04 

4.7546*03  -4.2896*04 

5.0006*02 

915 

544 

943 

696 

5.5236*03 

5.7526*03  -2.431E*04 

5.0006*02 

1453 

420 

694 

943 

-1.074E*04 

2.403E+03  -2.4946*04 

5.0006*02 

1442 

544 

714 

943 

1.5446*04 

6.8396*03  -5.3016*04 

5.0006*02 

1471 

420 

943 

714 

*1 .4536*02 

5.4906*03  -5.3246*04 

5.0006*02 

440 

474 

961 

449 

2.0A6E*03 

1.1576*04  -7.534E*04 

5.0006*02 

1460 

544 

449 

981 

9.050E+03 

8.270E+03  -4,744E*04 

5.000E*02 

1469 

474 

714 

981 

•1.3426*03 

1.3176*04  -4.5766*04 

5.0006*02 

1496 

544 

961 

714 

6.4826*03 

9.6456*03  -5.705E+04 

5.0006*02 

954 

565 

476 

449 

-3.874E*03 

1.720E+03  -6.8716*04 

3.0006*02 

1507 

474 

449 

476 

5.9406*03 

4.9486*03  -7.351E*04 

5.0006*02 

1516 

565 

514 

478 

*1.7426*03 

4.0916*03  -1.1206*05 

5.000E+02 

1525 

474 

476 

516 

7.5756*03 

1.1336*04  -9.7146*04 

5.0006*02 

972 

393 

1070 

467 

1.101E*04 

1.1126*04  -1.2916*05 

3.0006*02 

1534 

565 

467 

1070 

7.505E*03 

4.5946*03  -1.0426*05 

3.000E*02 

1543 

393 

514 

1070 

•4.663E*03 

5.534E*03  -1.398E+05 

5.0006*02 

1552 

565 

1070 

516 

•7.4506*03 

1.000E*03  -1.1436*05 

5.0006*02 
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CMD-POIHT  COm-LlMX  X  T  X,Y  VELOCITIES  X,Y  FORCES  MSS 


420 

109 

2.031E+00 

411 

96 

2.1326+00 

402 

120 

6.1S2E+00 

393 

227 

6.0936+00 

447 

0 

4.093+00 

429 

504 

4.135E+00 

456 

534 

2.0726+00 

465 

564 

6.1136+00 

474 

594 

4.043+00 

493 

0 

3.1036+00 

525 

0 

5.1176+00 

546 

0 

3. 0506+00 

585 

0 

5.0756+00 

576 

0 

4.105E+00 

642 

759 

3.131E+00 

615 

0 

3.076E+00 

495 

816 

2.099E+00 

633 

346 

5.143E+00 

433 

376 

6.131E+00 

651 

0 

5.093+00 

696 

933 

2.045E+00 

669 

0 

4.063E+00 

714 

990 

3.038E+00 

667 

1049 

6.092E+00 

516 

1079 

5. 064 E +00 

705 

0 

3.616E+00 

750 

0 

3.590E+00 

771 

0 

4.623E+00 

828 

0 

4.600E+00 

897 

0 

2.615E+00 

907 

0 

2.595E+00 

624 

0 

3.570E+00 

945 

0 

2.559E+00 

92' 

0 

5.633+00 

859 

0 

5.615E+00 

899 

0 

5.593E+00 

1031 

0 

4.593E+00 

963 

0 

2.539E+00 

981 

0 

3.550E+00 

678 

0 

4.562E+00 

1070 

0 

5.576E+00 

>STOP 

3.000E+00  -7.926E-02  4.314E-04 
6.992E+00  -1.894E-01  1.6466-01 
6.866E+00  2 >363-01  2.911E-01 
2.999E+00  3.754E-03  4.124E-02 
4.946E+00  3.223-01  9.7731-02 
6.926E+00  -2.463E-02  1. 393*01 
4.995E+00  1.003-01  -1.053-01 
4.9116+00  1.633-02  3.367E-02 
2.999E+00  1.473-01  4.260E-02 
5.954E+00  5.316E-02  -2.257E-01 
5.904E+00  1.350E-01  2.093-01 
3.9786+00  2.160E-01  -3.2051-02 
3.960E+00  3.623-01  4.203-02 
5.923+00  1.913-02  6.950E-02 
6.950E+00  -2.439E-01  2.264E  02 
4.963+00  1.097E-01  -1. 1591-01 
5.979E+00  -1.35E-01  -L.906E-01 
O.679E+00  5.093E-02  2.005E-01 
5.993+00  -6.8V4E-02  6.173-02 
4.931E+00  1.661E-01  2.523-02 
3.993+00  1.104E-01  -9.899E-02 
3.973E+00  1.043E-01  -7.978E-02 
2.999E+00  1.139E-01  1.451E-02 
3.953+00  3.827E-01  -1.253E-01 
2.WE+00  -1.103-02  -  3.503-02 
6.441E+00  -4.513-03  -7.033-02 
5.443+00  1.190E-01  2.097E-O1 
6.403+00  -1.665E-01  -4.274E-02 
5.429E+00  4.420E-02  1.521E-01 
6.465E+00  1.595E-01  1.291E-02 
5.470E+00  -1.313-02  2.526E-01 
4.463+00  1.743-01  -1.393-02 
4.490E+00  2.163-01  -6.608E-02 
6.391E+00  2.529E-02  -2.023-01 
5.405E+00  1.357E-01  6.610E-O2 
4.436E+00  1.349E-01  3.767E-01 
4.451E+00  -3.233-02  1.5886-01 
3.490E+00  1.458E-01  1.104E-01 
3.499E+00  1.379E-01  5.039E-01 
3.492E+00  9.105E-02  1.550E-01 
3.475E+00  -5.496E-02  1.040E-01 


5.391B+02  -4.995E+03  3.333E+02 
1.317E+03  -4.744E+03  3.333S+02 
-5. 9606+03  1.9036*04  3.3336+02 
2.343+02  9.801E+02  3.3336+02 
-1.4926+04  4.9396+03  1.3338+03 
5.513+03  6.349E+03  6.667E+02 
-8.405E+03  9.633+03  6,663+02 
-1.5166+04  4.9636+03  6.6676+02 
5.8116+03  1.7496+03  6.6676+02 
-1.1816+04  2.2126+04  1.3336*03 
6.919E+02  1.699E+04  1.3338+03 
4.3866+03  1.2876+04  1.3336+03 
1.564E+04  6.3666+03  1.333E+03 
2.133+03  3.3001+04  1.3336+03 
-3. 91 48+03  6.564E+03  6.663+02 
-5.304E+03  1.524E+04  1.3336+03 
2.903E+03  1.631E+04  6.663+02 
4.551E+03  2.713+04  6.663+02 
-5.930E+03  9.543E+03  6.663+02 
6.198E+03  9.5096+03  1.2336+03 
-3.749E+02  7.3456+03  6.663+02 
5.552E+03  -2.553E+03  1.3336+03 
-7.430E+02  2.059E+04  6.663+02 
-5.301E+03  -1.230E+04  6.667E+02 
2.562E+03  1.1356+04  6.663+02 
6.721E+03  -1.3Z7E+04  6.663+02 
1.352E+04  1.794E+03  6.663+02 
•4.783+03  1.5906+03  6.663+02 
2.903+03  -3,6916+03  6.663+02 
-8.959E+02  3.712E+03  6.663+02 
1. 6096+04  2.077E+02  6.663+02 
-1.304E+04  1.4206+04  6.667E+02 
1.1316+04  -4.133+03  6.663+02 
-9.123E+02  9.597E+03  6.667E+02 
-5.226E+03  8.272E+03  6.663+02 
-3.125E+03  2.541E+04  6.667E+02 
•2.7326+04  1.4566+04  6.663+02 
1.3126+04  1.6916+04  6.667E+02 
-4.6186+03  1.021E+02  6.663+02 
•9.6376+03  4.0936+03  6.667E+02 
9.4406+03  2.305E+04  6.663+02 
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Sample  Problem  No.  4 


A  projectile  hits  a  beam  and  breaks  It  Into  two  (fracture  based  on 
Griffith  theory). 


Figure  A4.  <jdec  5ampie  Problem  No.  4 
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START 

PROP  HAM  DENS 5 2000  K=i,QES  G-'j.OEfl  KNsi.E0fl  KS-l.EOA  F=0.1 
PROP  HAM  TF=3240, 

GRAVITY  0,  -10. 

ROUND  0.1 
DATP  .5  16.  NASS 
DA*'  .5  16.  INTERNAL 
FRAC  0.10 

6L0CK  ?1AM  C0NS=4  l.,0,  i .  >3 . 3b  5.  >3.35  j<  .0. 

SPLIT  0,1  6,1 
SPLIT  0,1.9  6,1.9 
SPLIT  1.5, -0,5  1.5,;,  5 
SPLIT  4. 5, -0.5  4. 5, 1.5 
DELETE  1.5, 4, ;>  0.1 
SPLIT  2. 5.1. 7  3.25.4 
SPLIT  3>2:  .1 ,6  4,4.1 
SPLIT  2. *5, l .27  3.75,2*4.3 
SPLIT  3.15,1.6  2.55.2.35 
Sr  LIT  2.6,3.')  4,3,05 
DELETE  i  ,2.5  1.9.3.J 
DELETE  3.75,5  1.9, 3. 3 
DELETE  2,o, 3. 4  1,9, 2. 4 
DELETE  2,3.75  3,15,3.3: 

ChANGE  1 ,5  1,2  SDEF 
LINK  3,3  4.1.5 
FIX  i,5  CM 

LOAD  1,5  2,3.35  -0-3E4  -C,95Eh 
VIB4  0  700  700  1400 
iiIND  0  6  0  7 


PLOT 

CYCLE  1000 

VIEU  700  1400  700  i-00 
PLOT 

CYCLE  200 

VIEW  0  700  0  700 

PLOT 

CYCLE  600 

VIEW  700  1400  0  700 

PLOT 

P  &  C  COP  DL  L 

STOP 

END 


PROBLEM  NO.  4  INPUT  FILE 
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>START 


PR0CRW1  UDEC:  UNIVERSAL  DISTINCT  ELB1ENT  CODE  (VERSION  1.2) 


THIS  IS  A  START  RUN 

>PR0P  MT<1  DENS=2000  M.0E8  C*1,0E8  KH=l*EOd  KS*1.E08  F=0.1 
>PRCP  MM  TF=3240. 

>CRAWTY  0.  -10. 

mie  o.i 

>DMf>  .5  16.  MSS 

STIFFNESS-DAHPINC  TERM  SET  TO  ZB® 

>D (tf  .5  16.  INTERNAL 
NOTE  -  ONLY  MSS-DOTINC  USED 
>FRAC  0.10 

JBLOC*  MM  CCNS=4  l.,0.  l.,3.35  5., 3.35  5,,0. 

>SPLIT0,1  6,1 
)SPLIT  0,1.9  6,1.9 
> SPLIT  1.5, -0.5  1.5, 1.5 
> SPLIT  4.5, -0.5  4.5, 1.5 
>DELETE  1.5, 4. 5  0,1 
> SPLIT  2.5, 1.7  3.25,4 
1SPLIT  3.25,1.6  4,4.1 
> SPLIT  2.65,1.87  3.75,2.48 
>SPLIT  3.15,1.8  2.55,2.85 
) SPLIT  2.6, 3.5  4,3.05 
>DEL£TE  1,2.5  1.9, 3.3 

DELETE  3.75,5  1.9, 3.3 
>DELETE  2.6, 3.4  1,9, 2.4 
>DELETE  3,3.75  3.15,3.35 
>CHANCE  1,5  1,2  SDEF 
)L1NX  3,3  4,1.5 
>FIX  1,5  0,1 

)*  SET  FORCE  OF  PROJECTILE 
>LOAD  1,5  2,3.35  -0.3E4  -0.95E4 
>VIEW  0  700  700  1400 
)UIND  0  6  0  7 
>PLOT 

) CYCLE  1000 

INITIAL  TINESTEP  '  6.325E-04 
CURRENT  CYCLE  COUNT  =  1000 
HI1EU  700  1400  700  1400 
>PLOT 

>CYCLE  200 

INITIAL  TIHESTEP  =  6.325E-04 
CURRENT  CYCLE  COUNT  =  1200 
>VIEU  0  700  0  700 
>PLOT 

>CYCLE  800 

INITIAL  TIHESTEP  *  6.325E-04 
CURRENT  CYCLE  COUNT  =  2000 
>VIEH  700  1400  0  700 
>P10T 

)P  6  C  COR  DL  L 
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MCI  OAM 


UQ  MT  C06t  CDflftMO  COORDS.  NASS 
1373  1  4  4.0Z3Et00  1.407EtO0  3.400Et03 

(SDET) 

011  012  021  022  Sill 

1.2MI*03  1.W4I-C3-1. 743-40  3.473E-03-4.3WEt03 

500  1  4  4.75OEtO0  3.000E-01  1.0CO£tO3 

(RIGID) 

193  1  4  3.140Et00  2.373EtOO  l.444£t03 

(RIGID) 

76  1  4  1.973Et00  l.403Et00  3.400Et03 

(8DGF) 

011  012  021  ED22  Sill 

W244E-03  1.WAE-03-1.7&3E-03  3.473£-03-4.3W£tC3 

l  1  4  1.250E*00  5.00CE-01  l-000£^33 

(RICID) 


WL.  non.  X,T, THETA  VEUX1TIES  X,t,I»ClA  FORCES 

1 . 443E*03  4.064E-02  4.201E-O2  1.545E-01  7.Z9&A3  2.U6Et04  l,135Et04 

SI  12  SI21  SIZ2  SA11  SA12  SA21  SA22 

2.yl5Et02-2.Vl5Et02-4.0V3£t80-3.025Et04-4.542£*0?-l.777Et03-2.231Et04 

1.M2EW2  O.OOOEtOO  O.OOOEtOO  O.OOOEtOO  4.203Et03  -2.357Et04  l.472Et03 

1.937E.02  -1.444E-02  -2.353E-01  -9.981E-03  l.7ilEt03  4.*43£t03  1.0S««02 

1.443E+03  -4.774E  02  -4.434E-02  -1.487E-01  -8.535Et03  2.070*04  -1.22^04 

SI  12  SI21  S122  SA11  SA12  SA21  SA22 

-2.T15Et02-2.Vl jEtO2-A.Ox3£tC3-3.025EtO4*4.542EtO2-l.777EtO3-2.Z3lEt04 

1.042E*O2  O.OOOEtOO  O.OOOEtOO  O.OOOEtOO  -4.721Et03  -2.4V2E*04  -1.5»4Et03 


P 


CONTACT  DATA 
CONTACT  NAT  COST 


X,T  COCMS. 

3.0MEt00  1 . 740E+00 
3.0O7EtOO  1.440Et00 
2.V44EtOO  1.7«5EtOO 
4.593EtOO  V.TV4E-01 
1.409EtOO  T.TV5E-01 
4.8V2£tOO  1 .OllEtOO 
l.llOEtOO  l.OlOEtOO 


FORCE 

NJChAL  SHEAR 
2.334E+03  2.334Et02 
1.067EtO4  1.047Et03 
5.255tt03  5.255E.02 
2.4lv£t04  -2.4l'tEt03 
2.322£t04  2.S22£t03 
O.OOOEtOO  O.OOOEtOO 
O.OOOEtOO  O.OOOEtOO 


DISPLACEMENT 
NORMAL  SHEAa 
-5.3171-05  -A.424E-04 
5.0T4E-05  -1.421E-02 
i.»02E-G4  -7.474E-03 
-2.374E-X  5.412E-02 
-2.423E-04  -4.533E-02 
2.1V5E-02  O.OOOEtOO 
2.400E-02  O.OOOEtOO 


F 


COWES  COORDINATES  *!N  X,Y  ORDES) 
RUO  1373 


3.0O2EtO0 

1.742£tOO 

4,»44EtOO 

1.T29EAX) 

5.053E+00 

1.034EtO0 

3,OJOEtOO 

4.4J1E-01 

BLOCK  500 
4.5001*00 

(.OOOEtOO 

5. OOOEtOO 

1. OOOEtOO 

5. OOOEtOO 

O.OOOEtOO 

4.SOOEtOO 

O.OOOEtOO 

BLOCK  195 

3.45ZItOO 

2.741EtOO 

3.373EtO0 

1.T23E.OO 

I'.HSEtOO 

I.724EtOO 

2.72OEtO0 

2.114EtOO 

BLOCK  74 
2.921EtOO. 

4.M5I-01 

♦.342E-01 

1.041E*» 

1.0l5EtO0 

l.iCVEtOO 

S.OlBEtOO 

1.743EtOO 

81(0  1 
I.SOOEtOO 

O.OOOEtOO 

1. OOOEtOO 

O.OOOEtOO 

:. OOOEtOO 

1. OOOEtOO 

I.SOOEtOO 

1. OOOEtOO 

DOMAIN  LINKED  LISTS 

DOMAIN  1189 

CONTACT 

1187,  BL0CX1 

1373,  BL0CK2 

195 

CORNER/EDGE 

CORNER 

544,  BLOCK 

1373,  R-LINX 

1393 

CORNER 

239,  BLOCK 

1373,  R-LIHK 

544 

CORNER 

106,  BLOCK 

1373,  R-LINK 

239 

CONTACT 

188,  BL0CK1 

1373,  BL0CX2 

500 

EDGE/CORNER 

CORNER 

120,  BLOCK 

500,  R-LINK 

558 

CORNER 

57,  BLOCK 

500,  R-LIW 

120 

CORKER 

532,  BLOCK 

500,  R-LINK 

37 

CONTACT 

846,  BL0CX1 

500,  BL0CK2 

1373 

CORNER/EDGE 

CORNER 

1393,  BLOCK 

1373,  R-LIW 

108 

CONTACT 

819,  BL0CK1 

76,  BL0CX2 

1373 

CORNER/EDGE 

CORNER 

709,  BLOC! 

76,  R-LINK 

215 

CORNS! 

1034,  BLOCK 

76,  R-LINK 

709 

CONTACT 

433,  BL0CK1 

1,  BL0CX2 

78 

CORNER/EKE 

CORNER 

334,  BLOCK 

1,  R-LIW 

96 

CORNER 

358,  BLOCK 

1,  R-LINX 

334 

CORNER 

21,  BLOCK 

1,  R-LINK 

358 

CONTACT 

144,  BL0CK1 

l,  6LOCX2 

78 

CORNER/EKE 

CORNER 

132,  BLOCK 

76,  R-LINK 

1034 

CORNER 

215,  BLOCK 

76,  R-LINX 

132 

CONTACT 

1429,  BL0CK1 

195,  BL0CK2 

76 

CORNER/CORNER 

CORNER 

1218,  BLOCK 

195,  R-LINK 

628 

CORNER 

1240,  BLOCK 

195,  R-LINK 

1216 

CORNER 

370,  BLOCK 

195,  R-LINX 

1240 

CORNER 

1405,  BLOCK 

195,  R-LINK 

370 

CORNER 

828,  BLOCK 

195,  R-LINX 

1405 

DOMAIN  937 (OUTER  BOUNDARY) 

CONTACT 

1187,  BLOCK 1 

1373,  BL0CK2 

195 

CORNER/EDCE 

CONTACT 

1425',  PL0CK1 

195,  BL0CT2 

76 

CORNER/CORNER 

CONTACT 

619,  BLOCK 1 

76,  BL0CK2 

1373 

CORNER/EDGE 

!  DOMAIN  870 

CONTACT 

848,  PL0CK1 

500,  8L0CX2 

1373 

CORNER/EDGE 

CORNER 

558,  BLOCK 

500,  R-LINK 

532 

CONTACT 

188,  BL0CK1 

1373,  BL0CX2 

500 

EDGE/CORNER 

1  DOMAIN  188 

CONTACT 

433,  BL0CK1 

1,  BL0CK2 

76 

CORNER/EDGE 

CONTACT 

144,  BLOCK 1 

1,  BLOCK! 

76 

CORNER/EDGE 

CORNER 

96,  BLOCK 

1,  R-LINK 

21 

i 

A36 

BLOCK  LIKED  LISTS 


BLOCK  1373 


CORBCR 

1393;  BLOCK 

1373,  R-U« 

106 

CONTACT 

419,  BLOCK  1 

76,  BL0CK2 

1373 

CORHER/EKE 

CONTACT 

1167;  H0CX1 

1373,  BL0CX2 

195 

CORNER/DCE 

cormr 

544;  BLOCK 

1373,  R-LDC 

1393 

cam 

239;  BLOCK 

1373,  R-LIME 

544 

cam 

106,  BLOCK 

1373,  R-LINX 

239 

CONTACT 

166;  BLOCK 1 

1373,  BL0CX2 

500 

EDCE/CORNER 

CONTACT 

646.  BLOCK! 

500,  BLOCK 

1373 

corner/edge 

BUCK  500 

COMR 

532,  BLOCK 

500,  k-UNX 

57 

CONTACT 

646,  BL0CX1 

500,  BLOCK 

1373 

CORtfR/EKE 

CORNER 

556,  HOCK 

500,  R-LINX 

532 

comer 

166,  BL0CX1 

1373,  BLOCK 

500 

EDCE/CORNER 

CORNER 

120,  BLOCK 

500,  R-LINX 

556 

CORNER 

57,  BLOCK 

500,  R-LINX 

120 

BLOCK  195 

cone* 

370,  BLOCK 

195,  R-LINK 

1240 

CORNER 

1405,  BLOCK 

195,  R-LINK 

370 

CORNER 

626,  BLOCK 

195,  R-LIMC 

1405 

CONTACT 

1167,  BLOCK 1 

1373,  BL0CK2 

195 

CORNER/EDCE 

CONTACT 

1429,  BL0CX1 

195,  BLOCK 

76 

CORNER/CORNER 

CORNER 

1216,  BLOCK 

195,  R-LINK 

628 

CORNER 

1240,  BLOCK 

195,  R-LINK 

1216 

BLOCK 

76 

CORICR 

709,  BLOCK 

76,  R-LMf 

215 

CORNER 

1034,  BLOCK 

76,  R-LINK 

709 

CONTACT 

433,  BL0CK1 

1,  BLOCK 

76 

CORNER /EDGE 

CONTACT 

144,  BL0CK1 

1,  BLOCK 

76 

C0RNER/EDC2 

CORNER 

132,  BLOCK 

76,  R-LINK 

1034 

CORICR 

215,  BLOCK 

76,  R-UNK 

132 

CONTACT 

1429,  BL0CK1 

195,  BLOCK 

76 

CORNEK/CORtCR 

CONTACT 

619,  BLOCKS 

76,  BLOCK 

1373 

CORNER/E DCE 

M.ont 

1 

CORtflR 

334,  BLOCK 

1,  R-LINK 

96 

CORNER 

356,  BLOCK 

1,  R-LINK 

334 

COOGR 

21,  BLOCK 

1,  R-LIW 

356 

CONTACT 

144,  BL0CX1 

1,  BLOCK 

76 

CORNEX/EDCE 

CORHER 

96,  BLOCK 

1,  R-LINK 

21 

CONTACT 

>5W 

433,  BLOCK 1 

1,  BLOCK 

76 

CORNER/EDGE 

rvt 


APPENDIX  B:  THREE-DIMENSIONAL  DISTINCT  ELEMENT  TEST-BED  CODE  (VERSION  1.0) 

USER'S  MANUAL 


Introduction 


1.  This  manual  describes  the  test-bed  code,  D3,  written  to  evaluate 
features  developed  in  the  design  of  a  new  three-dimensional  distinct 
element  program.  03  is  in  skeleton  form  with  several  facilities  provided 
for  in  the  code  but  not  completed  at  present.  The  input  commands  and 
program  operation  fol low  closely  those  given  for  UDEC. 


Input  Commands 


Notes:  Upper-case  letters  in  a  conmand  or  parameter  must  be  typed; 

the  remaining  letters  are  optional.  Lower-case  parameters  stand 
for  numeric  values.  Integers  must  be  given  for  parameters  start¬ 
ing  with  i,  j,  k,  1,  m,  n.  Real  numbers  may  be  given  as  integers, 
but  not  vice  versa. 

Input  is  free-format:  parameters  may  be  separated  by  any  number 
of  the  following  characters,  in  addition  to  spaces: 

»  (  )  .  / 

An  END  command  is  required  at  the  end  of  the  Input  fi’e  (after 
tne  STOP  conmand).  The  first  command  must  be  START  or  RESTART. 

*  *  comment  line 
+  *  continuation  line 


Block  Material  n  Constitutive  m  xl  yl  zl  x2  y2  z2  ... 

Create  a  rigid  block  of  material  number  n  and  constitutive  number 
m.  Defaults  are  n=l,  m=l,  if  m,  n  are  omitted.  The  block's 
surface  is  divided  into  triangular  faces.  Vertex  coordinates, 
(xl,  yl,  zl),  (x2,  y2,  z2),  etc.,  are  entered  three  at  a  time  for 
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each  triangular  face.  Continuation  lines  may  be  used  but  a  set 
of  three  vertices  defining  a  face  must  not  be  separated. 

Vertices  must  be  ordered  counterclockwise  looking  along  the 
outward  normal. 

CHange  xl  x2  yl  y2  zl  z2  Sdef  Material  n  Constitutive  m 

All  blocks  with  centroids  lying  within  the  range  xl<x<x2, 
yl<y<y2  and  zl<z<z2  are  changed  to  slmply-deformable  (Sdef)  or 
may  have  material  and  constitutive  numbers  changed. 

Cycle  n 

Oo  n  time-steps  (cycle  0  Is  permitted  as  a  check  on  data). 

DAmping  fcrit  freq  Mass 

Stiffness 

Internal 

Viscous  damping  is  applied  in  the  form  of  Rayleigh  damping, 
fcrit  Is  the  fraction  of  critical  damping  and  freq  Is  the  center 
frequency.  If  a  qualifier  is  not  given  as  the  third  parameter, 
full  damping  is  used.  The  word  "Mass"  eliminates  the  stiffness- 
proportional  dashpots,  and  "Stiffness"  eliminates  the  mass- 
proportional  dashpots.  The  word  "Internal"  causes  the  specific 
damping  to  be  applied  to  the  3  Internal  degrees  of  freedom  of 
slmply-deformable  blocks. 

Dump  n  m 

Dump  memory  to  printer  from  the  main  array  from  address  n  to 
address  m.  Internal  pointers  MFREE,  JUNK,  IBPNT  and  ICPNT  are 
also  printed.  MFREE  gives  the  highest  memory  location  that  Is 
currently  free. 

End 

Last  input  command. 

FRAction  f 

f  is  taken  as  the  fraction  of  critical  time-step  to  be  used. 
Gravity  gx  gy  gz 

Gravitational  accelerations  are  set  for  the  x-,  y-  and  z- 
directions. 

Print  Blocks  Faces  Velocities  VERtices  Contacts 

Data  are  printed  on  blocks,  faces,  block  velocities,  vertices  and 
contacts,  respectively. 
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PROperty  Material  n  keyword  value 
n 

The  first  parameter  must  be  the  specification  of  the  material 
number.  Material  properties  are  defined  for  material  number  n. 
Property  keywords  are: 


Bulk(orK) 

bulk  modulus 

G 

shear  modulus 

Density 

density 

KN 

contact  normal  stiffness 

K$ 

contact  shear  stiffness 

Cohesion 

contact  cohesion 

Friction 

contact  friction  coefficient 

JKN 

joint  normal  stiffness 

JKS 

joint  shear  stiffness 

Jcoh 

joint  cohesion 

JFric 

joint  friction  coefficient 

Restart 

The  program  Is  restarted  using  data  from  the  restart  file. 

RSet  v  ia  loff 

The  real  value  v  is  Inserted  in  the  main  array  at  address  la,  with 
offset  ioff. 

SAve 

The  current  problem  state  Is  saved  on  the  restart  file, 

STArt 

The  program  does  a  cold  start. 

Stop 

The  run  stops. 
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Prograa  Guide 


Parameters  and  Data  Group 


Offsets  for  block  data  array 


Hot*]  The  first  integer  in  «*ch  block  array 
—  (off tot  0)  it  the  block  typ*  nunter,  at  follow: 
t  rigid  block 
2  siaplydeforaable  block 

D  Pointer  to  next  block  in  block  1  ist* 

KF  Pointer  to  one  face  in  block's  face  list* 

KMT  Material  mater. 

ICONS  Constitutive  nuater. 

KBCOD  Code  number: 

0  free  block 
I  fixed  block 

KOI  Start  of  triple  pointer  to  x,y,z  coordinates 
of  block  centroid. 

KXD  Start  of  triple  pointer  to  x,y,z  coaponents 
of  velocity. 

ITD  Start  of  triple  pointer  to  x,y ,2  coaponents 

of  angular  velocity  (counterclockaise  positive) 
KV01  Block  volua*. 

KBH  Block  aass. 

XBI  Start  of  triple  pointer  to  aoaent  of  inertia 

about  x,y  and  z  axes. 

IffiFX  Start  of  triple  pointer  to  x,y,z  components 
of  block  centroid  force  sua. 

IBFT  Start  of  triple  pointer  to  x,y,z  coaponents 
of  block  centroid  aoaent  sua. 

KXL  Start  of  triple  pointer  to  x,y,z  coaponents 
of  load  applied  to  block  centroid. 

BO  Extension  pointer  (to  SDEF  data) 

IV  Pointer  to  on*  vertex  in  block's  vertex  list 
KC  Pointer  to  block's  contact  list 

Offsets  for  fact  data  array 


Mote:  The  first  integer  (offset  0)  contain* 
—  the  value  KFAC  to  denote  a  face. 


Off 

Pointer  to  next  face  on  this  block. 

M 

Pointer  to  host  block. 

r+rrt 

m 

Pointer  to  first  connecting  face. 

KF2 

Pointer  to  second  connecting  face. 

1 


IP3  Pointer  to  third  correcting  fact# 

IW1  Pointer  to  first  vertn  of  this  face. 

W2  Pointer  to  oecond  vertex  of  this  fxr. 

«V3  Pointer  to  third  vertn  of  this  face. 

Offsets  for  vertn  data  array 

Motet  The  fit  tt  integer  (of fret  0)  contain* 

* —  the  value  MR  to  denote  a  vertn. 

BW  Pointer  to  next  vertn  on  thie  block 
IVX  Start  of  triple  pointer  to  a^z  coordinates 
of  vertn. 


Muter  of  writ  in  dm  arrays 


MAL  Block 

WFR  Face 

MMt  Vertex 

NVCH  Contact 

IM)C  Degenerate  contact 

Array  Unit* 


HTOP  Sin  of  Min  array  (IA). 

MMT  Raxiaua  number  of  aaterials. 

NCONS  Raxiaua  constitutive  nunbers. 

HTTP  limber  of  block  types  (rigid ,  SDEF,  etc.) 

Head  codes  (contents  of  first  integer  in  data  groups) 


MIC  -  1  Rigid  block 

HSDEF  3  2  Sieply-defomable  block 

HFAC  Face 

IWER  Vertex 

HCON  Contact 

Contact  codes 


HCEE  edege-to-edge 
HCVF  vertex -to- face 
HCDC  degenerate 

Rain  Canon  Black  Variables 


LtlC(BO)  Buffer  for  current  input  line  in  Ai  forati 
LINEKB6)  Buffer  for  next  input  line. 

IPHT(I)  Pointer  to  start  of  paraeeter  1  in  LDS(  ) 

after  renova  I  of  blarits;  etc. 

(WHO)  vector  of  zero  length 

ERFLAC  .TRUE,  if  an  error  has  a  cured. 

STFLAC  .TRUE,  if  the  first  input  line  has  been  processed. 
COFLAC  .TRUE,  if  the  current  line  is  a  contiuation. 

HCFLAC  .TRUE,  if  the  next  line  is  a  continuation. 

JRPSAV  Index  of  last  conputed  GOTO  in  HON. 

IGW  Error  nuater. 

AM  Pointer  to  list  of  spare  aeeory  poupt. 

IVREE  First  unused  aeeory  address. 

I  BLOCK  Current  block  nuater. 

1ST AC*  Stack  pointer. 
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tore 

Currently  requested  nutter  of  eye  In* 

NCTOT 

Toni  nutter  of  cycln* 

7ML 

Tine-step. 

FRAC 

Requested  fraction  of  critical  tine-step. 

I  ROl/TE 

Routing  nutter,  usad  in  ain  routine. 

HOC 

Output  lint  count* 

WAGE 

Output  page  count* 

JMPC&i 

Routing  nutter  for  continuation  lino  in  GEM* 

ALPHA 

test  duping  coefficient* 

BETA 

Stiffness  duping  coefficient. 

CONI 

Batping  factor  <1.0-AlPHteTUL/2.0> 

C0N2 

Batping  factor  <1.0/<1.0+ALPm*TMlj?.0» 

m 

BETA/TD0. 

an 

Internal  eats  datpinq  coefficient  for 
sieply-defornable  blocks. 

ClB 

Duping  factor  <1.0-ALPB»TDEL/2.0) 

C2& 

Duping  factor  <1.0/0. OMlPteTDEL^.O)) 

DBCRAD 

PI/180 

PI 

3.14159 

ATOL 

Distance  tetieen  particles  at  teich  a  real  contact 
is  forted. 

cm 

Dm 

Distance  between  particles  at  Uich  a 
degenerate  contact  it  forted 

Em 

Fm 

IBPHT 

Pointer  to  list  of  blocks. 

law 

Pointer  to  list  of  contacts. 

AKN(I) 

Mortal  contact  stiffness,  oaterial  I. 

AKS(I) 

Shear  contact  stiffness,  uteris)  I. 

MU<1) 

Contact  friction  coefficient,  nterial  I. 

GOH(I) 

Contact  cohesion,  Nterial  I. 

AKNJ(I) 

Joint  norul  stiffness,  Nterial  I. 

AKSJ(I) 

Joint  shear  stiffness,  Nterial  1. 

AHJJ(I) 

Joint  friiction  coefficient,  Nterial  I. 

COHJ(I) 

Joint  cohesion,  Nterial  I. 

DOS(I) 

Density,  Nterial  I. 

BULK  I ) 

Bulk  eodulus,  Nterial  I. 

S>CMt<I> 

Sheer  eodulus,  Nterial  1. 

ALAMl(I) 

Lane  constant,  Nterial  I. 

ALAM2(I> 

Lane  constant,  Nterial  1. 

CWV<3> 

vector  of  pavitational  acceleration. 

A(  ) 

tein  array. 
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Two  tetrahedral  blocks  are  created.  The  lower  block  Is  fixed, 
and  the  upper  block  allowed  to  come  Into  contact  with  the  upper  vertex 
of  the  fixed  block.  Gravity  acts  In  the  -2  direction.  Since  the 
centroid  of  the  upper  block  Is  not  directly  above  the  fixed  vertex,  the 
block  translates  and  rotates,  and  develops  shear  forces  at  the  contact 
as  well  as  a  normal  force. 
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a.  initial  State 


HW»  NAM  0*2000  DM  IGA  IS».5E0  F*1.0 
BUCK  NAM  (0,0,0)  (1,0,0)  (1,1,0) 

♦  (1,0,0)  (1,0,1)  (1,1,0) 

♦  (0,0,0)  (1,1,0)  (1,0,1) 

♦  (1,0,0)  (0,0,0)  (1,0,1) 

BLOCK  MM  (0, -0.5,1)  (1.5, -O.S,l)  (1.5, 1,1) 

♦  <1.5, -0.5,1)  (1.5, *0.3,3)  (1.5, 1,1) 

♦  (0, -0.5,1)  (1.5, 1,1)  (1.5, -0.5, 3) 

♦  (0, -0.5,1)  (1.5, *0.5, 3)  (1.5, -0.5,1) 

FIX  0  1  0  1  0  l 

CRAW  0  0  -10 
CTC  1 

PRINT  CON  VEL 
CYC  5 

PRINT  CON  m 
CYC  10 

PRINT  CON  VEL 

STOP 

EM) 


INPUT  FILE 


program  #3i  tot-kb  s-o  iutdct  uwr  com 


W1S  IS  A  STMT  MM 

IfUDP  MM  #’2000  DMEI  1X.5E8  M.O 
>NOCX  MM  (0,0,0)  (1,0,0)  (1,1,0) 

>♦  U,0,0>  (1,0,1)  (1,1,0) 

>♦  (0 fifi)  (1,1,0)  (1,0,1) 

>♦  <1,0,0>  (0,0,0>  (1,0,1) 

>IU)CI  MT»1  (0, -0,5,1)  (1.5, -0.5,1)  (1.5, 1,1) 
>♦  (1.5, -0.5,1)  (1.5, -0.5, 3)  (1.5, 1,1) 

>♦  (0, -0.5,1)  (1.5, 1,1)  (l.S,-0.5,3> 

>♦  (0, -0.5,1)  (1.5, -0.5, 3)  (1.5, -0.5,1) 

>FU  0  I  0  1  0  1 
>«AM  «  0  -10 

iac  i 

TMESTEP  =  2.542E-0* 

awfflff  cycle  cant «  i 
>mw  cm  va 

contact  type  stocx-i  wjxi-2  x<i 


104  (vstra-mr: 

177  (DEEEMERATc) 

MjOCK  code  uxmjt 


09  1.000£*00  O.OOOE*00  1.000E.O0  O.OOOE.OO  0.000£*00  0.0005*00  4.4475*01 


E  UXMJT  UrtCT  UZDOT  TOOT  TYOOT  TZWT 

0  O.OOOE+OO  0.000£*00  -2.5CZE-03  0.0005*00  0.0005*00  0.000£*00 

1  O.OOOE*CO  0.0005*00  0.0005.00  0.0005*00  0.0005*00  0.0005*00 


0.0835-00  -0.0035-00  -1.5095-02  4.9791-05  4.9795-05  0.0005*00 
0.0005*00  0.0005*00  0.0005*00  0.0005*00  0.0005*00  0.0006*00 


TIICSTEP  •  2.5025-04 
CURRENT  CYCLE  COHT  *  4 

>FRIHT  CW  VEL 

CONTACT  TYPE  EL0CI-1  NOa-2  X(l)  X(2)  ((3)  FS(1> 

104  (VERTEX-EDCE)  1  09  l.OOOE^O  0.0005*00  1.0005*00  4.424-01 

177  (DEGENERATE)  69  1 

sua  coos  uxmjt  irroor  uzbot  too:  tybot 

09  0  0.0635-00  -0.0035-00  -1.5095-02  4.9795-05  4.9795-05 

1  1  0.0005*00  0.0005*00  0.0005*00  0.0005*00  0.0006*00 

>CTC  10 

TDMSTEP  •  2. 5025 -U4 
CURRENT  CYCLE  COUNT  •  14 

JPRMTCCN  Via 

CONTACT  TYPE  IUtt-1  UOB-2  X(l>  X(2)  X 

104  (VBTTEX-CK5)  1  09  1.0005*00  0.0005*00  1.000C* 

177  (0— tUK)  09  1 

tun  con  ukdot  unor  uzbot  toot  tywt 

09  0  1.3075-05  -1.3075-05  -3.3951-02  9. 1415-04  9.1415-04 

1  1  0.0005*00  0.0005*00  9.0005*40  O.WNiOO  9.0005*00 


9X2) 

-4.4245-01 


no)  m 

0.0005*00  1.3791*03 


1.0005*00  9.0001*00  1.0005*00  2,4535*01 


71(2) 

•2.4535*01 


n<3>  m 
0.0005*00  0.1495*03 


ssesissesaEr 


